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Abstract: We present an extension to next-to-leading order in the strong coupling 
constant g of the AMY effective kinetic approach to the energy loss of high momen¬ 
tum particles in the quark-gluon plasma. At leading order, the transport of jet-like 
particles is determined by elastic scattering with the thermal constituents, and by in¬ 
elastic collinear splittings induced by the medium. We reorganize this description into 
collinear splittings, high-momentum-transfer scatterings, drag and diffusion, and par¬ 
ticle conversions (momentum-preserving identity-changing processes). We show that 
this reorganized description remains valid to NLO in g, and compute the appropriate 
modihcations of the drag, diffusion, particle conversion, and inelastic splitting coeffi¬ 
cients. In addition, a new kinematic regime opens at NLO for wider-angle collinear 
bremsstrahhmg. These semi-collinear emissions smoothly interpolate between the lead¬ 
ing order high-momentum-transfer scatterings and collinear splittings. To organize the 
calculation, we introduce a set of Wilson line operators on the light-cone which deter¬ 
mine the diffusion and identity changing coefficients, and we show how to evaluate these 
operators at NLO. 
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1 Introduction 

Jets are a key observable in the relativistic heavy-ion program [1-4]. Advances in re¬ 
constructing jets at the LHC [5, 6] challenge our ability to understand the difference 
between jet development in the hot medium created in a heavy ion collision, compared 
to development in the vacuum or near-vacuum environment of a proton-proton colli¬ 
sion. While early theoretical studies concentrated on understanding the leading hadron 
in a jet (see ref. [7] for an overview), the more inclusive jet reconstructions which are 
now possible experimentally demand a theoretical description of the full jet evolution, 
including the evolution of all radiated daughters. 

Several groups have put forward modeling frameworks for doing this [8-12]. It is 
fair to say that these approaches have some commonalities. Generally they separate the 
excitations into high-energy partons associated with the jet, and low-energy partons or 
a scattering medium with a characteristic energy scale T (the local medium tempera¬ 
ture). Then, one attempts to follow the evolution of the high energy partons, which 
will eventually create the hadrons reconstructed as a jet. The jet partons are considered 
to interact with the medium in two important ways. They scatter elastically, and they 
are induced to radiate or split. Different frameworks differ in whether both possibilities 
are considered, and in exactly how the splitting processes are computed (how is long¬ 
distance coherence handled? Is the radiated daughter assumed to have a small fraction 
of the energy? What model for the medium interactions, and what other approximations 
are made?). 

Typically the division of processes into distinct types - here elastic scattering and 
inelastic radiation - is justified at leading order, but at subleading orders they often 
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cannot be clearly distinguished. What happens to the treatment of jet-medium in¬ 
teraction at subleading order? Is it possible to pursue a next-to-leading order (NLO) 
calculation, in the sense that the elastic and splitting interactions between the jet par- 
tons and the medium are treated beyond leading perturbative order? ^ In this paper 
we explore this question by extending a framework where it is clearly posed - the 
AMY/McGill/MARTINI approach [10, 20-24], where there is a clear power-counting 
prescription for determining what is leading and subleading order for the jet-medium 
interaction. The approach starts with the assumptions that the medium is thermal and 
weakly coupled, so the interactions between jet partons and the medium can be com¬ 
puted in thermal perturbation theory. One also assumes that the medium is thick, such 
that the formation times of processes under consideration are shorter than the scale of 
variation of the medium. This approximation has sometimes been criticized, and it can 
be improved upon without overturning the rest of the approach [25]. But for the pro¬ 
cesses which will be most interesting here - processes involving small momentum transfer 
or intermediate opening angles - the scattering or formation times are relatively short, 
so this should be considered a separate issue. 

The philosophy of the framework is as follows. We follow one or more “hard” 
approximately on shell partons traversing the medium. We assume that the medium 
has a local temperature T 3> Aqcdj and distinguish a parton as hard if its energy 
E satisfies exp(—£'/T) <C 1. Particles failing this criterion are assumed to join the 
thermal medium; but no attempt is made to track the back-reaction on the medium 
properties [26]. 

The jet parton evolution and jet-medium interactions are dictated by finite tem¬ 
perature perturbation theory. Whereas vacuum perturbation theory is an expansion in 
the strong fine-structure constant ag = jAir, this expansion is spoiled by soft-particle 
statistical functions Usioj rsj gT) r-u 1/g entering in thermal Feynman graphs. These soft 
contributions must be resummed to obtain a hnite leading order answer, and give rise 
to subleading corrections suppressed by a single power of g. We have recently shown 
how to compute these subleading corrections in the context of hard real [27] and vir¬ 
tual [28] photon production. Here we extend that treatment to the case of jet-medium 
interactions. 

As a scattering environment, the essential attribute of QCD (or any gauge theory) 
is that there is a large cross section for small-momentum-transfer scattering processes. 
These are responsible for the high rate of particle splitting. They also cause complica¬ 
tions when including elastic scattering, since they give a large rate of small momentum 

^Some authors have used the term NLO in a different sense: for instance that the initial parton pro¬ 
ducing processes or the final fragmentation processes are treated at NLO, though the medium interactions 
are still leading order [13, 14], or within the Higher Twist formalism [15, 16[, or that higher-order, double- 
logarithmic corrections to the jet-quenching parameter are considered and resummed [17-19]. Here by 
NLO we mean a beyond-leading-order treatment of the way the jet interacts with the medium. 
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exchanges, both in the transverse and longitudinal components of the momentum. At 
next-to-leading (NLO) order, new processes arise, which can be understood physically as 
overlap and interference between sequential scattering processes and as scatterings with 
the emission or absorption of soft {E ~ gT) excitations [29]. These contribute both to 
transverse momentum broadening and to longitudinal momentum loss and broadening. 
They are most easily computed in a way which does not cleanly separate them into elas¬ 
tic and inelastic processes, and indeed it is not clear that the distinction is important 
or well posed. And they overlap with the infrared limits of both the elastic scattering 
and the splitting processes. However, frequent and small momentum exchanges need 
not be separately identified and tracked. In traversing enough medium to significantly 
modify a jet, the jet partons will undergo several such soft processes, in which case 
a statistical description should be sufficient. This motivates an approach in which we 
give a Fokker-Planck (Langevin) description of soft scatterings, as drag and diffusion 
processes. 

The philosophy of our approach will therefore be the following. We will introduce 
infrared scales /r_L, /r||. All scattering and emission processes which change a jet parton’s 
momentum by more than /x_|_ [| will be handled explicitly. All processes which change 
momentum by less than this scale, including the NLO effects alluded to above, will be 
incorporated as momentum diffusion and drag coefficients, which can be neatly defined 
in field theory as correlators of field strength operators on light-like Wilson lines. We will 
compute these drag and diffusion coefficients at the NLO level, as well as providing an 
NLO accurate procedure for computing the larger-transfer elastic and splitting processes. 
We will show explicitly how to perform a matching so that the choice of the scale g drops 
out in the final results. 

The drag and diffusion coefficients account for momentum exchange with the medium 
through soft gauge-boson exchange. Soft fermion exchange with the medium can change 
the identity of a quark to a gluon and vice versa. We call such identity changing pro¬ 
cesses conversion processes, and introduce a medium coefficient (analogous to the trans¬ 
verse momentum broadening coefficient q or the energy loss e) which parameterizes this 
conversion rate. As with the drag and diffusion coefficients, all identity changing scatter¬ 
ing processes with momentum exchange greater than g will be treated explicitly, while 
identity changing scatterings with small momentum transfer are incorporated into the 
conversion rate. The conversion rate will be defined as a correlator of soft fermionic 
operators on light-like Wilson lines. 

The calculation of the processes involving soft momentum transfers (drag, diffusion, 
and conversions) and of the corresponding light-cone Wilson line correlators requires a 
resummation scheme known as the Hard Thermal Loop (HTL) effective theory [30, 31], 
which is the QCD analog of the Vlasov equations [32]. These formalisms are well known 
to be computationally complex, and at first sight, any calculation beyond leading order 
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in the coupling would seem extremely challenging (see [29] for an example). However, 
Caron-Huot has shown [33] that HTL correlators (and statistical correlators more gen¬ 
erally) simplify greatly when computed at light-like separations, which are exactly those 
that must be evaluated to determine the energy loss, diffusion, conversion, and collinear 
radiation rates of highly energetic particles propagating in a plasma. 

Intuitively, these simplifications can be seen to arise because the energetic partons 
are propagating almost exactly along the light cone. Hence they are probing an essen¬ 
tially undisturbed plasma, at least as far as the soft, classical background is concerned. 
Informally, we can say that this background “can’t keep up” with the hard particles 
traversing the plasma. Thus, the soft correlations that the latter probe are statistical in 
nature rather than dynamical. Those simplifications, as we shall show, are at the base 
of the NLO extension being presented. 

A pedagogical review of these recent developments in the understanding of HTLs 
has been presented by two of us in [34]. There the main results of this paper, i.e. the 
reorganization of the kinetic theory we have mentioned, as well as the results of the 
computations to NLO of the needed rates and coefficients, have been partly anticipated. 
Due to the review nature of [34], the presentation there has been more pedagogical and 
most details and technical aspects have been omitted for the sake of brevity and clarity. 
Here we will present the detailed derivation of the reorganization of the kinetic approach, 
as well as the explicit calculations of the coefficients and rates. Furthermore, [34] was 
limited to a plasma of gluons only, again for ease of illustration. We advise readers unfa¬ 
miliar with Hard Thermal Loops, and especially with the recent developments discussed 
before, to explore [34] first. 

The paper is organized as follows: in Sec. 2 we present the LO framework in the 
standard formulation, which divides into elastic (2 -f-)- 2) and inelastic (1 -H- 2) processes. 
Readers familiar with that approach can skip directly to Sec. 3, where we introduce 
our reorganization in terms of large-angle scatterings, diffusion, conversion and collinear 
processes. In Sec. 4 we give an overview of the NLO corrections, which are dealt with 
in detail in Sec. 5 for collinear processes. Sec. 6 for diffusion. Sec. 7 for conversion and 
finally Sec. 8 for the semi-collinear processes, which first contribute at NLO and smoothly 
interpolate between the other three. A summary is presented in Sec. 9, together with 
our conclusions. Extensive technical details are to be found in the appendices, such as 
the NLO calculation of longitudinal momentum diffusion. 

The paper is rather long and detailed. Some readers will be more interested in 
learning how to apply its results to their own (numerical) treatment of jet modification, 
without necessarily following all the details. In the online version of the paper, we have 
put boxes around the key equations which must be included in an implementation of 
our results. These readers can skip most of the text and focus on these boxed equations. 
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2 The leading-order kinetic approach 

Our aim is to track the time evolution of a small number of highly-energetic jet-like 
particles as they propagate through a medium. We will refer to energies and momenta 
of order E as hard, of order T as thermal and of order gT as soft^. The hard particles are 
very close to the mass shell, with energy ^ p ^ E and virtuality \po — p| < g^T. We 
will assume that this energy is large enough that exp{—E/T) 1 and can be neglected, 
but we will not treat T/E as an explicit expansion parameter. Thus, for instance, we do 
not distinguish between a rate that is of order g^T and one that is g'^y/TE. Moreover, we 
will often find convenience in using light-cone coordinates, specifically those defined by 
the hard four-vector P. If, without loss of generality in an isotropic medium, p points in 
the £ direction, then for a generic vector K we can define k~ = kP — and k~^ = . 

This normalization, already adopted in [27], is nonstandard, but we find it convenient 
because dk^dk^ = dk~^dk~, and because we will frequently encounter cases in which 
k~ ~ 0, in which case k^ k, k^ k, k'^ with our conventions. The transverse coordinates 
are written as k_L, with modulus k±. 

Let us start from the effective kinetic theory developed in [21]. The Boltzmann 
equation reads 

(^ + V • V,) r{p,x,t) = -cr'[/] - (2.1) 

where f°‘{p,x,t) is the phase space distribution for a single color and helicity state 
quasiparticle of type a (/“ = dN°‘/{dPxdPp)). In the collision operator, at leading order 
in the coupling g, one needs to account for 2 o 2 and effective 10-2 processes. The 
2 o 2 rates are given by the simple 2 0-2 diagrams of QCD, such as those shown in 
Fig. 1, which also establishes our graphical conventions. The 1 o 2 rates describe the 
collinear radiation from the jet-like particles, which is induced by multiple soft scatterings 
with the background plasma, see Fig. 2. Although apparently suppressed by powers of 
g, multiple scatterings contribute at leading order under the provision that: (a) the 
momenta of the hard lines are nearly on shell and collinear to each other (i.e. 0 < g, 
where 9 is the emission angle^), and (b) the momenta K of the soft gluons are space-like 
k'^,kx ~ gT and k~ < g^T. A complete leading order treatment of collinear radiation 
must consistently resum these soft scatterings to account for the Landau-Pomeranchuk- 
Migdal (LPM) effect [20, 35-38]. 

^ The notation and terminology used here is summarized in App. A, and closely follows our previous 
work [27]. 

® In the case where P and Q are both thermal, such as when dealing with the thermal photon rate, 
then the angle is of order g. In the case of interest, i.e. P hard, there are two different possibilities. If 
either Q or P — Q are thermal, i. e. there is a hierarchical separation between the emitted particles, then 
the angle is again of order g. If instead the splitting is more democratic, with no hierarchical separation, 
then the angle can become as small as gpTjE. 
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Figure 1. Typical diagrams contributing to 2 ^ 2 processes at LO. Double lines represent 
hard or thermal particles, which have at least one momentum component of the order of the 
temperature or larger. Parallel double lines without arrows can be either gluons or quarks. 
When particle identities need to be specified, quarks are identified by the fermion flow arrow 
and gluons by the curly line. In all diagrams in the paper, time is understood to flow from left 
to right. 



Figure 2. A typical diagram contributing to 1 -O- 2 processes at LO. The single curly line is a 
soft gluon. The crosses represent the soft thermal scattering centers - see, for instance, ref. [20] . 


In detail, the collision operator reads (dropping for brevity the spacetime depen¬ 
dence, which is local) 


cr^[/](p) 


and 
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{27r)'^ + K - P' 


x{r{p)f\k) [l±r{p')] [l±Ak')] 

-np')f\k')[i±rip)] [i±f\k)]], 


K') 


( 2 . 2 ) 
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Y1 / '^p ~p ~^bc{p\p'p, q'p) 

2\p\^l^a^Jo 

X {rip) [i±f\p'p)] [i±riq'p)] - f{p'p)nQ'p) [i±/“(p)]} 

+ ^(IpI+ 7ab(p'p;p>9p) 

\P\ Jo 

X [rip) f\qp)[i±rip'p)] - fip'p) [i±rip)][i±f\qp)]}, m 

where the sum runs over the species 6 c((i) in the scattering/splitting event, and the 
splitting kernel 7 ^^ is defined in Eq. (5.1-5.4) of Ref. [21], see also Eq. (5.1). We are also 
using the shorthand notation 


for the Lorentz-invariant integration. The matrix elements M and transverse-momentum 
integrated matrix elements 7 will be discussed in the following. = 2dR is the degen¬ 
eracy of the particle o: two spin degrees of freedom and (Jr color degrees of freedom, 
where Jr is the dimension of the representation of a. For quarks it is dir = Nc, for 
gluons “ 1 - 

The hard particles are very dilute, and therefore we only need to track the inter¬ 
actions of these modes with the thermal and soft constituents. This can be done by 
defining 5f 

rip, X, t) = n“(p, r(a:, t), u(®, t)) + (5/“(p, x, t), (2.5) 

and linearizing the Boltzmann equation in this quantity. Here n is the (local) equilibrium 
distribution, written generally as a function of the local temperature T and flow velocity 
u. In the following we will work in the local rest frame where n becomes the Fermi-Dirac 
distribution n^(p) or the Bose-Einstein distribution nsip)- Substituting Eq. (2.5) in the 
collision operator and dropping terms which are of order yields 


= 


J\p\T^a 


Y, [ Mf^ip,k-,p',k')^ i27r)U^^\P + K - r - K') 

bed 


X |d/“(p) n^(/c) [l±n'^(p')±n'^(/c')] — d/'^(p') n'^(A:') [l±n^(A;)] 

-n'^ip')6fik')[l±n\k)]}, ( 2 . 6 ) 


and 
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c^^^WKp) 


X] / ^(\P\~P' ~ ^bc{P'^ P'P^ Q'p) 

2\p\^Ua^J, 

X {'5/“(p) [l±n'’(p')±ri'=(g')] - W\p'p)n''iq') + n\p')5f{k'p)]^ 

+ ri^5 ilbip'mp^qp) 

\P\ Jo 

X {^/“(p) n'(g) - 6np'p)[l±n\q)]}. (2.7) 

These are to be used in a linearized Boltzmann equation for the hard components 5/“ (p ), 

• V,) 5r(p,rr,t) = -CLO[^[/] = -Cl^^[6f] - Cr'[5/]. (2.8) 

In the 2 o 2 collision integrals, soft glnon and fermion exchanges must be screened to 
avoid logarithmic divergences. At leading order, the bare t— and tt— channel propagators 
may be replaced with their Hard Thermal Loop counterparts in the infrared to render 
the collision integrals hnite. This procedure (which is detailed in Appendix A of [39]) 
provides the leading weak-coupling description for soft exchanges, and is correct to order 
for hard exchanges. However, while this regularization prescription provides the 
correct leading order answer, it is not easily generalized to NLO. Further, the approach 
mixes different physics at different scales. In the next section we will re-examine the 

2 -f-)- 2 collision rates, incorporating soft t, u—channel exchanges into drag, diffusion, and 
conversion coefficients, which cleanly reflect the physics of the Debye sector. Then, in 
Sec. 4, we will compute these transport parameters at NLO. 

3 A reorganization of leading order: large-angle scattering, drag and 
diffusion, and conversions 

The leading order picture we have just described, with distinct 1 2 collinear processes 

and 2 •(-)• 2 processes dressed with HTLs for IR hniteness, starts to be ill-dehned at NLO. 
Consider the soft limit of the 2 •(-)• 2 processes. In the case of a soft gluon exchange, as 
shown in Fig. 3, we obtain a process which changes the hard four-momentum P by a 
small amount Q ~ gT, without changing the particle identity. We call snch a process 
a diffusion process, since, as described in Sec. 3.2, they can be treated in a diffusion 
approximation. 

In the case of a soft quark exchange, as shown in Fig. 4, one obtains an identity¬ 
changing process. Here a hard gluon with momentum P is turned into a quark with 
an almost equivalent momentum, up to 0{gT). We then call these processes conversion 
processes, and we will deal with them in a different way, inspired by the NLO thermal 
photon rate [27]. 


- 8 - 





Figure 3. The soft limit of a t— or u—channel gluon exchange diagram. P is the hard momentum 
and Q is the soft gluon momentum. 


uQ 


Figure 4. The soft limit of a t— or m— channel quark exchange diagram. P is the hard momentum 
and Q is the soft quark momentum. 


Now consider a collinear 1 -H- 2 process in the limit where one of the hard/thermal 
legs becomes soft^, as shown in Fig. 5. In the first graph, the soft gluon emission 



Figure 5. The soft-iti limits of a 1 O 2 process. The diagram on the left amounts to a diffusion 
process at NLO, whereas the diagram on the right amounts to a conversion process. 

contributes to the (longitudinal) diffusion of the hard particle. Similarly the soft quark 
emission contributes to the hard quark conversion rate. At NLO we will then need 
to subtract these limits from the collinear 1 •(-)• 2 region and treat them as part of the 
diffusion or conversion processes respectively. 

'^LPM interference is suppressed in this case [27]. 
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To summarize, at leading order we can rewrite the right-hand side of Eq. (2.8) as 


= - 

Here is the 2 0 2 collision operator restricted to large momentum transfers, Q S> 

gT. Dehning this scattering rate requires regularization procedure, which we describe 
in the next section, Sec. 3.1. notates a diffusion approximation to the collision 

integral for small momentum transfer, Q ~ gT. This is discussed in Sec. 3.2, where 
the LO longitudinal and transverse diffusion coefficients are extracted from the screened 
2 -H- 2 rates. Similarly, notates the conversion processes, and the appropriate LO 

conversion coefficients are found in Sec. 3.3. The precise value of these diffusion and 
conversion coefficients depends on the regulator, but the dependence on the regulator 
cancels to leading order when C'large^ ^difr^ ^conv taken together in Eq. (3.1). 
Finally, (7®°*^ consists of the collinear 1 •(-)• 2 rates after excluding (or subtracting) 

the diffusion and conversion-like emissions shown in Fig. 5. These soft emissions (which 
were originally included in the rates) are limited in phase space to K gT, 

and their exclusion constitutes an 0{g) correction. Thus, at leading order = 

C^‘^‘^[df]. We therefore will present the explicit form of only when we describe its 
NLO corrections in Sec. 5. 

3.1 Large-angle scattering 

In this section, we describe the integration of the 2 -f-)- 2 matrix elements with large 
momentum transfer, Q 3> gT, which enters in the leading order collision kernel 
in Eq. (3.1). This is completely straightforward, but integrals of the bare matrix ele¬ 
ments must be regulated with some scheme. The cutoff regulator chosen in this section 
conveniently matches with the calculations of the diffusion and conversion coefficients 
in Sec. 3.2 and 3.3. 

In more detail, to evaluate one needs integrate the matrix elements listed 

in Table 1, i.e. the standard, leading-order QCD matrix elements, summed over all color 
and spin indices, with the Mandelstam variables s = —{P + K)‘^, t = —{P — P')‘^ and 
u = —{P — K')"^. To regulate gluon and fermion exchanges in the t and u channels we 
use the integration technology of Ref. [39] , which treats each channel differently. 

Singly-underlined matrix elements come from gluon exchange diagrams, and are 
those that, in the soft limit, give rise to gluonic IR divergences, corresponding to diffusion 
processes. Similarly, doubly-under lined matrix elements come from fermion-exchange 
diagrams and give rise, in the same limit, to conversion processes. To illustrate the 
regularization scheme, let us consider the contribution from the scattering of different 
quark species qiq 2 ^ qiq 2 , which is given by the square of a single t-channel diagram. 


- cr^'W]- 


(3.1) 
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Table 1. Squares of vacuum matrix elements for 2 -o- 2 particle processes in QCD-like theories, 
summed over all spins and colors. qi and q 2 represent fermions of distinct flavors, qi and q 2 are 
the associated antifermions, and g represents a gluon. 


The < 71^2 •t-t contribution to reads^® 

^4 /•+00 p2p—uj 




g 


/ -\-oo pzp—LJ r 

doj dq 

-oo Jo J (c 


dq I dk6{q — |w|) 

{q-u})/2 


(27r)3 IQp^ J_ 

XI {p)np{k)[l — np{p — u) — np{k + u)] 

— [df^^ {p — uj) np{k + w) + 5/'^^ {k + w) np{p — u 


Jo 27r 


)] [l-n,(k)|},(3.2) 


where the techniques of [39] have been followed, by (i) eliminating one of the three 
integration variables in Eq. (2.6) with the momentum-conserving (5-function, (ii) shifting 
one of the remaining ones to q = p — p' = k' — k, (hi) introducing uj = p — p' = k' — k, 

® N.b. the sum over c and d in Eq. (2.2) yields a factor of two. 

® For ease of illustration, we are considering for large-angle scatterings a simplified case where 5f is 
a function of p rather than of p. Details on the phase space integration in the latter case can be found 
for instance in [40]. 
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and (iv) performing the angular integrations. The remaining angle (j) represents the 
azimuthal angle between the {p,q) plane and the {k,q) plane. 

The Mandelstam variables become 


s 



(P + p'){k + k') + q^ - cos{(j))y/{App' + t){Akk' + t) , 


t = u}'^ — (3.3) 


and = p', k'^ = k' imply 


UJ-pq = 


co^-q^ 


2p 


u — k ■ q = 



(3.4) 


It is then easy to see how the unscreened logarithmic divergences show up for w, g ~ 
gT <C k,p. To separate off the divergent region (which will match with the diffusion 
operator described in Sec. 3.2), we change integration variables from u!,q to 

w, gj_ with g_i_ = = \Jq^ — We can then place an IR cutoff T 3> Pq^ S> gT on 

g_L, leaving 



Eq. (3.5) implicitly depends on the cutoff pq^. For small q± the dominant u region is 
also small, and the fermion distribution can be approximated as, np{p — a;) ~ np{p). 

In the small q± regime, the scattering rate integrated over u turns out to take a very 
simple form in terms of this variable, which is the real motivation for its use. In addition, 
the physical interpretation of q±_ in this regime is the transverse momentum transferred 
to the p particle; specihcally, in terms of the P-dehned light-cone coordinates, we have 
for soft Q 


~2 2 
Q± = Q± 


q 


+ 


1 + ^ + 0 


m) 




ql fg^T^ 
= a;-k — -hO ' 


2p 


p^ 


(3.6) 


Finally, let us analyze the power counting. Above the cutoff, when angles are large, 
the contribution to the collision operator is of order g^T, up to powers of T/E. When 
q,uj ^ gT, the ())-averaged matrix element is proportional to p'^k'^/q'^ ~ up to cor¬ 
rections, which combined with dqduj ~ g^ and with another g"^ coming from the expansion 
of the curly brackets for small Q make the singly-underlined exchanges contribute to LO 
in the soft region, with a In(gf) enhancement. The same happens (without cancellations) 
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for the doubly-underlined matrix elements. Non-underlined matrix elements with st or 
tu at the denominator are suppressed by a further power of in the soft region. Since 
the integration is finite, can be pushed to zero there for simplicity. Singly underlined 
matrix elements with a v? at the denominator present the same divergences; they can 
be dealt with by swapping the k' and p' labels and using the same parameterization. 
Matrix elements with at the denominator are not sensitive to the soft region; hence, 
at leading order, they can be integrated without cutoffs as well. 

Fermion exchanges, and in particular the log-divergent doubly-underlined t- or u- 
channel exchanges, can be treated with the same techniques and cutoffs. For illus¬ 
tration, the t-channel quark exchange contribution to qiqi gg scattering is 


(3.7) 


The cancellations of the leading IR behavior in the gluon exchanges, as well as the 
matching to the diffusion and conversion processes will be dealt with in the next sections 
and in App. D. 



3.2 Diffusion processes 

In this section we will describe the diffusion collision kernel, in greater detail. The 

cumulative effect of a large number of small momentum-transfer collisions that preserve 
the identity of the hard particles can be summarized by a Fokker-Plank equation [41, 42] 


cf^W] = 


A 

dp'’ 


VD{p)p"^r{p) 


iA!_ 

2 dp'dpi 


p'fqLip) + 


fp>)q{p) <5/“(p) 


(3.8) 

App. D directly shows how the diffusion operator arises at leading order from the 
screened 2 •(-)• 2 collisions kernel, Eq. (2.2). There are three coefficients that enter in 
this effective description: q is the standard transverse momentum broadening, is the 
longitudinal momentum broadening and is the drag coefficient. They are dehned as^ 


^d(p) = -^^, g(p) = ^((Apx)^), ffi{p) = j^{{/\pLf), (3.9) 

^These coefficients depend on the species a. However, as we shall show, to leading and next-to-leading 
orders in g this dependency reduces to a simple Casimir scaling in the representation of the source a, so 
we drop this label in the text for simplicity. 
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where pi^ and p± are the longitudinal and transverse components relative to the large 
momentum p. 

These coefficients can be determined through the interaction rates [41-43], i.e. 


II 


(3.10) 

9r(p) = 


(3.11) 

q{p) = 

f ,2 2 c^r(p,p-hq) 

(3.12) 


where r(p,p± q) is the transition rate from initial hard momentum p to final hard 
momentum® with q soft. A regulator which cuts off the Q integrations is implicit, 

and the values of these coefficients will in general depend on the chosen scheme. Rather 
than determining r(p,p± q) and evaluating the integrals in these equations directly, it 
is convenient (especially at NLO) to use field-theoretical definitions for the coefficients 
in Eq. (3.8) 

The transverse scattering rate dT /d?q± at large momentum is traditionally param¬ 
eterized hy C{q±) 


lim 

p^oo 


dr{p,p + q±_) 
(Pq± 


C{q±) 

(27r)2 • 


(3.13) 


In p —)• oo limit the hard particle’s behavior eikonalizes, and C{q±) can be defined 
in terms of a specific Wilson loop [33, 44] in the plane (for propagation in the 

positive z direction). Using this Wilson loop definition, C{q±) and q have been evaluated 
at leading [45] and next-to-leading orders [33]. In particular, at leading order the result 

is 


CR{qx) = 9 ^Cr I ^^27rd(g° - q,)G--(Q) = q^CrT 


rUr 


qliql+ml) 


(3.14) 


where R labels the representation of the source and 
leading order Debye mass, q then reads at LO 


ml = g^T^{N,/3 + Nf/6) 


9 = y = 9^CrT 


d?q± rnl ^ g^CRTml 

(27r)2 (q2 ) 27r ^ rriD ’ 


(3.15) 


where, since we are in the p —)• oo limit, q_\_ = q± and we have used pq^ as UV regulator. 

What makes the Wilson loop definition particularly attractive is that it can be 
evaluated [33] using the (much simpler) Euclidean, dimensionally-reduced Electrostatic 
QCD (EQCD) [46-50]. This made the NLO computation possible [33], and opened the 

®Since the exchanged momentum is soft by construction, there is no ambiguity in the identification 
of the hard outgoing line 


- 14 - 















door to recent non-perturbative lattice measurements [51, 52]. These formal definitions, 
as well as those for related light-front operators, are summarized in App. B of [27] and 
reviewed in [34]. 

These techniques and results, like most eikonal expansions, are based on a large 
momentum expansion, p 3 > T or gT. In App. D we study the finite-p corrections, 
showing that T/p suppressed corrections are really corrections in gTjp] and the first 
correction involves vanishing odd integrands, so the first nonzero corrections from this 
expansion are 0{g‘^) even for p ^ T, and are therefore irrelevant at the level of precision 
we are seeking here. Therefore we can use the leading (and later, subleading) order 
calculations in the strict Wilson-line limit which we have just discussed. 

To fully specify the diffusion operator in Eq. (3.8) we also need to evaluate the 
longitudinal diffusion and drag coefficients, and rjo- To this end, we will first compute 
the diffusion coefficient and then use fluctuation-dissipation relations to determine 
the drag (see below). At the practical level, we introduce a Wilson-line based definition 
for qi^ in the p —)• oo limit, or equivalently at leading order in T/p.^ In App. C we will 
give a more formal justihcation for our definition, whereas in App. D we show that, as 
in the previous paragraph, hnite-momentum corrections start at 0{g‘^), and are thus 
irrelevant to current accuracy. 

Intuitively, longitudinal momentum diffusion occurs because the longitudinal force 
along the particle’s trajectory has a nonzero correlator. Experience with q and heavy 
quark diffusion [53, 54] suggests that q^ should be given by a lightlike longitudinal 
force-force correlator. The force is determined by the electric field in the direction 
of propagation, which motivates the following operator definition for in the large 

momentum limit: 

„2 ^- 1-00 

= — / dx+Tr ([/ij(-oo,0,0_L;x+,0,0_L)h^Ui,F^‘'(3:+,0,0_L) 

J — OQ 

C/ij(x+, 0, Ox; 0,0, 0 x)F^‘"(0)C/r(0; -oo, 0, Ox)hpU,) . (3.16) 

Here v = (1,0,0,1) and v = (1/2,0,0, —1/2) are null vectors that are chosen to maintain 
our light-cone conventions (he. p'*' = —v ■ P, p~ = —v ■ P), and F^’^Vfj^Vu = F~^ = 

is the electric field along the propagation direction. We are using a matrix notation, so 
that F^’^ = and Ur is a straight Wilson line in the representation R of the 

source ^ 

UR{Y-X) = Pe^p(-ig ds{Y-X)-A r{s{Y-X) + X^ . (3.17) 

Gauge fields and matrices in the Wilson lines are both to be understood as path ordered. 
Eq. (3.16) comes from the eikonal approximation, he. the replacement of the highly 

®In the following, q and Ql are understood to be in the infinite-momentum limit unless otherwise 
specified. 
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energetic particle with momentum p with a Wilson line in the appropriate representation 
along its classical trajectory. We also note that this definition of has the correct 
“amplitude times conjugate amplitude” structure required to enter in a rate. 

We now evaluate Eq. (3.16) at LO: we simply contract the two F fields, obtaining 
a forward Wightman correlator, i.e. the diagram shown in Fig. 6, which reads 



Figure 6. The leading-order soft contribution to q^. The Wilson lines before and after the two 
black dots, which represent the vertices, cancel at leading order, whereas the one between 
the two dots always turns into an adjoint line, which we have represented as a double line. The 
curly line is a soft HTL gluon. 

Ql = 9^Cr (3.18) 

where G^{Q) is the HTL-resummed forward propagator and the integral is understood 
to run over soft momenta only. The integration sets q~ to zero and, as we show 
in App. D, brings this expression into agreement with the one obtained from the rate- 
based dehnition in Eq. (3.11). Note that only the even-in-q"'" part of G^{q'^, q~ = 0, q±) 
contributes to the integral. Then, using the fluctuation-dissipation theorem, G^{Q) = 
(1 + nB{q^))p{Q) with p{Q) = Gr{Q) — Ga{Q), we expand for small q^ = q'^ ^ gT to 
hnd 

qL = g^GR j ^^^^;^y|^rg+(G^"(g+,g_L) - G;^"(g+,g_L)), (3.19) 

up to an 0{g‘^) correction. Numerical integration is straightforward, using the HTL 
propagators given in App. B. Beyond leading order, however, one would be plagued 
with intricate multi-dimensional numerical integrals. However, as we anticipated in the 
introduction, we can perform the q~^ integration (and similar ones elsewhere) by resorting 
to the analyticity sum rule techniques developed in [27, 33]. Since retarded (advanced) 
two-point functions are analytic in the upper (lower) half-plane in any time-like or light- 
like variable, we can deform the integration contours away from the real axis onto Cr 
(k"*"! gT, Img+ > 0) and Ca (k’*"] gT, Img+ < 0), as depicted in Fig. 7. Along 

The following derivation has been anticipated in [34]. 
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Im 



Re 


Figure 7. Integration contour in the complex integration, and the deformation we use to 
render < 7 + ^ gT. Gr runs above the real axis and Ga below. 

the arcs the longitudinal and transverse propagators simplify greatly, i.e. 



(3.20) 


where = 171^^12 is the gluon asymptotic thermal mass. The end result is then 



(3.21) 


where contributions smaller than l/q~^ in Eq. (3.20) are not needed, as they would 
only give rise to power-law terms in the cutoff on q'^ which would then cancel against 
contributions from larger scales. As in the q case, we have used as a transverse 
regulator, since the {q~^, q±) and (cu, q±) coordinates differ by 0{g‘^) (all 0{g) corrections 
vanish under integration). 

The sum rule we have just obtained is the bosonic equivalent of the one presented in 
[27]. Let us remark that the longitudinal and transverse contributions to G^~{Q) contain 
poles at q~^ = g“/2 ± iq± {q^ = 0), which, being on both sides of the complex plane, 
appear to violate analyticity. However their residue cancels in the sum of longitudinal 
and transverse components. As observed in [33], they are artifacts of the decomposition 
into Lorentz-variant longitudinal and transverse modes and their contribution has to 
vanish in all gauge-invariant quantities. 
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We also remark that the same result (3.21) has been obtained in a different way 
in [55] for energy loss, which is related by an Einstein relation. As shown there, once 
the difference in regularization between q± < and q < fig is taken into account, 
Eq. (3.21) agrees with the numerical results of Braaten and Thoma [43] for u —1. 

Having determined q^, the drag coefficient r]iy{p) is constrained by the requirements 
that the Fokker-Planck description be equivalent to the Boltzmann one and that interac¬ 
tions with the medium tend to drive the hard excitations towards equilibrium [42, 56, 57]. 
Since we have taken a classical particle approximation for the hard particles, the equi¬ 
librium form is 6f{p) oc exp(—p/T). The drag (in a given regularization scheme) is 
determined from qL(p) and q(p) by adjusting the value of Pnip) so that Eq. (3.8) ap¬ 
proaches equilibrium, i.e. its right-hand side vanishes for 6f{p) oc exp(—p/T). Since q 
and are p-independent up to 0{g‘^), the equilibration condition yields to following 
relation: 


Vd{p) 


Ql 

2Tp 


+ 


2p‘ 




(3.22) 


The consistency of this condition is verified by direct computation of rjoip) and q^ at 
leading order in App. D. Inserting this relation between the coefficients into the diffusion 
equation, Eq. (3.8), we find 




' Sfjp) 2T + pd6f{p) I d^dfip) - 

Tp 2pT dp^ 2 dip^Y 


Ql - 


1 d6f{p) 
2p dp^ 


+ -V^ 

^ 4 P-L 


^f{p) 


(3.23) 


which is our final form for the diffusion operator. 

We end by making a few remarks about the equilibrium condition and The first 
term in Eq. (3.22) comes from the simple Einstein relation that arises in the infinite 
momentum limit, i.e. q^ = —2Tdpi^/dt + 0{l/p). The relative 0{l/p) terms can then be 
obtained by imposing equilibration on Eq. (3.8). In App. D we will show how, at leading 
order, those 1/p terms can be determined explicitly, how the diffusion picture matches 
exactly with [5/] at large Q and how different cutoff schemes can be implemented. 
It is also worth stressing that the 1/p terms to Eq. (3.22) and equivalently to rfn do 
not come from a, T/p expansion, but only from the 5 <C 1 expansion. Up to relative 
there are no l/p'^ terms. Einally we remark that, in the simpler case where 6f is a 
function of p rather than p, as in footnote 6, the contribution proportional to q vanishes 
in Eq. (3.23). 
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3.3 Conversion processes 

The conversion-process part of the collision operator can be simplified as 


cr’Wi =sr-(p)Tfr,(p) - «/“(p)^r“4-(p), 


Nf 


/^conv 

^9 


i«/i=E 


2=1 




dp 

dA 





(3.24) 

(3.25) 


(3.26) 


representing a rate for each species to disappear due to conversion to another type, 
and a rate for that species to appear due to the conversion of another type to the 
type in question. The conversion rates describing these processes can depend on 

momentum p, and they also implicitly depend on the regularization scheme. They do 
not, however, depend on the exchanged momentum at leading or next-to-leading order. 
To see this, consider Eq. (3.7) for LO,q± ~ g. One has that the statistical factors, once 
expanded for (7 <C 1 , yield 


i^6f‘^^{p)np{k)[l+ nB{k)] - 5/®(p) ns(A;)[l - np.(A;)]| ^1-^0 . (3.27) 

Similarly, as we shall show in more detail in App. D.2, the HTL-resummed and (p- 
averaged matrix elements, once expanded for small Q, are to leading order even in cj, 
up to 0{io/T,uj/p) corrections.^^ Furthermore, the {oj^qp) and ((?^,(?_l) coordinates are 
equivalent up to another odd-in-w correction, as shown in Eq. (3.6). All these odd, 
subleading corrections vanish upon duj integration, so corrections first arise from uP‘ jT'^ 
type corrections or the product of two w/T corrections, which are both safely NNLO. 

At leading order the rates can simply be obtained from the aforementioned even- 
in-cu term in the HTL-resummed, doubly underlined matrix elements, whereas at next- 
to-leading order soft-gluon loop corrections need to be considered. To this end, we 
find it convenient to define the conversion rates in terms of gauge-invariant Wilson line 
operators, following the work on the soft contribution to the photon rate in [27], where 
we showed that the leading and next-to-leading order soft contributions were obtained 
from similar operators. Physically, the amplitude for a quark to convert to a gluon 
involves a quark propagating in from an early initial time, and being converted to a 
gluon by the insertion of a quark destruction operator ip. The rate is the product of this 

^^Non-underlined fermion exchange matrix elements, such as u/s in q\g -Pp qig scattering, are sup¬ 
pressed by two powers of g. 
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amplitude with its conjugate, and we must integrate over the time difference between 
the quark annihilation event in the amplitude and in its conjugate. Eikonalizing, the 
propagation of the quark turns into a fundamental Wilson line, while the gluon which 
propagates between the earlier and later ij}, ip insertion is represented by an adjoint line. 
This leads to the following Wilson-line representations for the conversion processes: 

o2 r+oo 

J dx+(w[[/i.(-oo,0,0x;x+,0,0x)r“^(x+,0,0u)^ 

X Ua{x+, 0,0; 0,0, 0_l)V’(0)T^C/f(0; -oo, 0, 0_l)] ) , (3.28) 

0^ f+°° 

^g°^qip) = J (Tr[C/^(-oo,0,0x;x+,0,0_L)T“V’(a;+,0,0_L)^ 

X 0,0; 0,0, 0±)'iIj{0)T'^Ua{0-, -oo, 0, 0_l)] ^ , (3.29) 

r^-°-(p) = r“-(p), r“-(p) = r“-(p). (3.30) 


The traces appearing in the rates are over the Dirac and color indices. At leading order 
the rates read 




g'^Cp f (fq± g^Cprul^ /ig^ 

Ap J (27r)2 gqj-p ’ 


dp 

dA 


■pconv 


ip) 


LO 


(3.31) 

(3.32) 


m 


2 = g^CpT'^lA 


IS 


where we have used the light-cone sum rule obtained in [27, 58]. 
the asymptotic mass of quarks. The regulator is the same used in the large angle 
and diffusion regions. In App. D.2 we show how the evaluation of the appropriate part 
of the HTL-resummed 2 •(-)• 2 collision operator in this momentum region leads to the 
same result. 


4 Next-to-leading order corrections: Overview 

The reorganization we have presented in the previous section allows us to introduce 0{g) 
corrections to the collision operator. For convenience we identify two different sources, 
i.e. loop corrections and mistreated regions. The former arise by adding a soft gluon 
loop to a diagram, which, in the finite-temperature power counting, gives rise to an 
0{g) contribution. The latter instead originate from integrating over 0{g) regions of 
the leading-order phase space where one particle becomes soft, without being treated 
correctly as an HTL quasiparticle. One such example is mentioned at the end of Sec. 3 
and in Fig. 5, where a soft, final-state gluon in a 1 2 process gives rise to a finite 


- 20 - 















contribution to the LO 1 •H' 2 collision operator. As we shall show, this indeed represents 
an 0{g) region of the 1 -H- 2 phase space; its evaluation, as well as the evaluation of all 
such mistreated regions, requires the identification of the limiting behavior of the LO 
calculation in that region. Such behavior will then have to be subtracted from the 
proper, HTL-resummed, evaluation of that region, which, in the example of Fig. 5, will 
be done when dealing with at NLO. 

In the large-angle region, loop corrections are suppressed by a factor of as long 
as the momentum transfer stays large. But the LO evaluation, in the form of Eqs. (3.5) 
and (3.7), mistreats the region where an incoming gluon is soft. This region will be 
properly addressed in the semi-collinear region, which we shall introduce later on. We 
defer other considerations on the necessary subtraction to that point and to Sec. 8.^^ 

In the collinear region, we will encounter both loop corrections and subtraction 
regions. The former arise from adding extra soft gluons to the scatterings that broaden 
the hard particles, inducing their splitting. They correspond to the NLO corrections 
to C(gj_) [33], which have been already mentioned after Eq. (3.13). The asymptotic 
masses of the hard particles also receive 0{g) corrections that contribute at NLO. In 
Sec. 5 we will discuss in detail those corrections, as well as three mistreated regions; the 
aforementioned overlap with the diffusion region, an altogether equivalent one with the 
conversion sector and finally one with the semi-collinear region. 

In the diffusion sector, Eq. (3.8) remains valid to NLO. Its coefficients go-, Ql and q 
all receive 0{g) loop corrections. Those to q are known [33]. In Sec. 6 we will set up the 
calculation of the 0{g) corrections to q^, through the field-theoretical definition (3.16) 
and the causality-based sum rules. The details of the evaluation will be presented in 
App. F. It requires the subtraction of a mistreated 0{g) region in its LO evaluation, as 
well as of the aforementioned diffusion limit of the collinear sector. Finally, go can be 
determined through the equilibration condition (3.22). 

In the conversion sector, the operators defined in Eqs. (3.28)-(3.30) receive 0{g) 
loop corrections from the addition of one extra soft gluon. In Sec. 7 we will show how 
these operators are equivalent up to NLO to their abelian counterparts. Hence, the 
0{g) corrections can be extracted from the soft-sector contribution to the NLO photon 
rate in [27]. In this case too there are subtractions from mistreated regions in the LO 
conversion and collinear rates. 

Einally, a new kinematical region enters at NLO, the aforementioned semi-collinear 
region. It corresponds to medium-induced splittings with larger virtuality, transverse 
momenta and opening angle, respectively of order gT'^, gT^ and ^/g^^ Other differences 
with respect to the collinear region are that the kinematics now allow the soft gluons to 

At the NLO level there is also a linear in fi± divergence in evaluating q, which is canceled by a linear 
in fj.1- soft-gluon effect in the hard scattering regime, see [33, 59]. This divergence and mistreatment 
simply cancel; so will not discuss it further, directing the interesting reader to those papers for details. 

^^Up to respective factors oiT/E and sjTjE in a democratic splitting case, similarly to footnote 3. 
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be either space-like or time-like (hence the overlap with the soft limit of the large-angle 
region), that LPM interference is suppressed and that the soft gluons can change the 
small minus component of the hard/thermal particles’ momentum. We will deal with 
this sector in detail in Sec. 8. As mentioned, we will have to subtract the mistreated 
overlap regions of the large-angle and collinear regions. 

We conclude this overview by sketching the form of the NLO corrections: 


6CaW] =6C2°^\6f] + 6C^^[6f] + [^Z] 


(4.1) 


Here and in what follows 6 refers to an NLO contribution. The first term consists of 
the loop correction to the collinear sector. In the second term the form of the diffu¬ 
sion equation (3.23) remains unchanged, but the parameters, q and ql, receive NLO 
corrections from soft loops. In particular, the corrections to the longitudinal diffusion 
coefficient, depends logarithmically on an ultraviolet cutoff, Similarly, 

the momentum dependence of the conversion rates remains unchanged, oc 1/p, but the 
overall magnitude of the rate depends logarithmically on This dependence on the 

ultraviolet cutoff in the diffusion and conversions collision kernel cancels in the complete 
kernel when the semi-collinear emission rates are included. In the semi-collinear case, 
^NLO ggpygg ^g Infrared cutoff limiting the semi-collinear emission of soft quarks and 
gluons. 

To compute each of the collision operators in 6C, the phase space regions which were 
mistreated at LO must be subtracted as counterterms. This replaces the mistreated LO 
terms with the full NLO result, and generally removes power divergences in soft loop 
integrals: 


5C^^[5f] =AC^^[5f] - - SC^LsnuAm, (4-2) 

5crAm =AcrAm - ^c^^cZisuuAm - sc^^An^sn^Am, a-a 

6c~^Am =Acr“^-“'V/] - scizAXtAm - (4.4) 

In each case, the subtraction terms arise from a mistreatment in a specific region of phase 
space from one of the four LO collision kernels in Eq. (3.1). For example, in the first 
line treats the diffusion process with NLO accuracy by including the appropriate 

soft loops, while the two counterterms arise because the LO collinear and LO diffusion 
collision kernels give incomplete contributions to the diffusion process at NLO. 

We will devote the next four sections to evaluating in turn the four contributions to 
the NLO collision operator given in Eq. (4.1). 


5 The collinear region 

Here we discuss the NLO corrections, and subtractions, needed to establish splitting 
processes to this order. But for completeness and context, and to set notation, we begin 
by presenting the leading-order result. 


22 - 



5.1 Leading-order recapitulation 

At LO which is [ 22 ] 


(5.1) 


(5.2) 


where the 0 -function multiplying the last term prevents a double counting of the gg final 
states (equivalently one may use a ^ symmetry factor). = T{p,puj) is the rate 

for a particle a with hard momentum p to emit (u > 0 ) or absorb (a; < 0 ) a gluon 
(quark in the case F^g) with energy (longitudinal momentum) 6fq^q is either <5/^. or 
6fqp. Eq. (5.1) applies both for quarks and antiquarks, provided a consistent labeling of 
LO in T^qq is chosen. At leading order these rates read [20, 22]^® 




where x = oo/p is the momentum fraction of the outgoing gluon or, in the qq final state, 
of one of the two fermions, h = p x g is the two-dimensional invariant describing the 
transverse separation of the final states. Note that for QCD, in the cases q ^ qg and 
g ^ qq Cr = Cp, whereas for g ^ gg Cr = Ca- F(h,p, w) determines the transverse 
evolution of the system; it is to be determined through an equation which resums multiple 
soft interactions. In momentum space it has the form of an integral equation, whereas in 

keeping with the notation in the other sections, we label lj the longitudinal component of one of 
the outgoing momenta. 

^®The distribution functions in [22] are summed over spin, color and flavor, so that the factors of dA/dp 
and dp/dA vanish in Eqs. (5.1) and (5.2). However, the g ^ qq rate in [22] and subsequent references 
(see [10, 24]) was missing the factor of dp/dA that appears in Eq. (5.3). Indeed, the group-theoretical 
factor for this process should read Cpdp/dA = Tp = 1/2. 
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position space it is a differential one. The latter will be described in App. E; the former 
reads [ 20 ] 


2h = i5E{h,p,co)F{h) + J 


d?k_i_ CF{k±) 
(27r)2 Cf 


(CH-CV2)[F(h) -F(h-wkx)] 


+^[F(h) - F(h + pkF)] + ^[F(h) - F(h - (p - a;)k^)] 


(5.4) 


For the case of 5 —)• qq, {Cr — CaI‘^) multiplies the term with F(h — pk_i_) rather than 
F(h — wkx). The equation depends on two inputs, CF{k±) and 6E{h.,p,uj). The former 
is C{k±) for a fundamental source, whereas 5E is the energy difference between the initial 
and final collinear particles. It reads 


6E{h,p,u}) 




+ 


m“ 


2pu}{p — cj) 2a; 


+ 


^oop—u) 
2{p - Uj) 


mlop 

2p ’ 


(5.5) 


where m^p is the asymptotic mass of the particle with momentum p, as summarized in 
Eq. (B.3). 


5.2 The collinear sector at next-to-leading order 

Subleading corrections to collinear splitting are treated, for the case of photon produc¬ 
tion, in [27]; the case here is conceptually similar. We must identify any NLO corrections 
to splitting for generic kinematics; and we must identify any limits of the kinematics 
which contribute an 0{g) faction of the total splitting rate, but which overlap with the 
kinematics in another region we are studying. 

The NLO corrections for generic kinematics enter as two corrections which arise 
when solving Eq. (5.4), specifically, 0{g) corrections to C{q±) and to the asymptotic 
masses entering in 6E. The computation of these masses to NLO has been carried out 
in [60] using Euclidean techniques and is reviewed in [27, 34]. It is only due to soft gluons 
and depends on the nature of the particle (quark or gluon) through a simple Casimir 
scaling: 

(5.6) 

The NLO collision kernel SC{q±) has been computed in [33], as a Hrst application of 
the mapping to the Euclidean theory. All one needs to do to treat generic momenta at 
NLO is to include these two corrections into Eq. (5.4). We review how to do so, using 
impact-parameter-space methods, in App. E. 


= -9^Cf 


Trrii 

2f 


5Ml = -g^CA 


Trrii: 

2f 


5.3 Subtraction regions 

Besides these generic-momentum corrections, there are also corners of the collinear- 
splitting kinematics where it starts to overlap with other processes - momentum diffu¬ 
sion, identity change, and 2-^2 scattering. Each regime represents an 0{g) suppressed 
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fraction of the total contribution from splitting processes, so a correct leading-order 
treatment is sufficient. Unfortunately, in each regime at least one approximation made 
in arriving at Eq. (5.3), Eq. (5.4), Eq. (3.14), or Eq. (5.5) breaks down. We handle this 
in two steps. First, we find out what contribution the (naive) leading-order splitting 
calculation actually contributes in each region. Then, we perform a more complete NLO 
calculation of the specific kinematic corner of interest, subtracting the (naive) leading- 
order splitting contribution we have found, since it is already incorporated via the LO 
splitting treatment. The remainder of this section carries out the calculation of the LO 
splitting behavior in each kinematical corner. 

In each relevant corner, 6E ^ f d^k±C(k±), so that, physically, the formation time 
l/SE of the collinear particles becomes much shorter than the time between collisions, 
estimated by (f d^k±C(k±))~^. In this case emission amplitudes associated with different 
scattering events become incoherent, and it is sufficient to treat emission as a sum of 
the rate arising from each scattering event (LPM suppression is small), up to 0(g) 
corrections which we can neglect. In the diffusion and conversion cases this happens 
because the denominators in Eq. (5.5) become smaller by a factor of g, whereas in the 
semi-collinear case becomes larger by 1/g. Therefore, we first obtain the generic 
solution for 6E S> f d^k±C(k±): if we solve Eq. (5.4) by substitution as in [27, 61] we 
have at leading order^® 

-2h 

ImF(h,p,a;) = —---, (5.7) 

dE[h,p, CO) 

which in turn yields 


ReF(h,p,a;) = 


/ 


d^k± CF{k±) f 


+- 


dE(h,p,uj) J (27r)2 
h h.+ p kx 


<5E(h) <5E(h+pkx) 


+ 


Cf 

Ca 

2 


1 


Cr- 


Ca \ h h — a; kx 

) [<5E(h) “ <5E(h-a;kx)_ 


h 


h — (p — u;)kj 


(5E(h) (5E(h — {p — a;)k_i_) 


(5.8) 


5.3.1 The diffusion limit 

Let us now specialize to the soft gluon region, which corresponds to the diffusion limit. 
Explicitly, one has oj —)• gT in the q gq and g gg processes. In the case of the latter 
process, there is also a (p — ca) ~ gT region which appears in the loss term in Eq. (5.2) 
but is absent from the gain term. Since its contribution is identical (the loss term and 
Eq. (5.3) are symmetric around ca = p/2) this compensates for the relative factor of two 
between g ^ gg gain and loss terms in Eq. (5.2), yielding for Eq. (5.2) a limit of the 

The g -F)- qq case, which has a different color structure in curly braces, is not dealt with explicitly. 
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(5.9) 


form of the ^i-proportional part of Eq. (3.23)^'^. We then have 

= 2^ + ^ ^ ’ 

which indeed is of order gT and larger than the collision operator by a factor of l/g. 
Eq. (5.8) can be integrated over cPh, symmetrized and expanded for small uj to become 


(f'h 


= Sp^Cax^{ 1 — 2x) 


2h • Re F(h,p, w) 
(fq± 


soft g 


( 2 vr)^ 


f d?kx Cpikx) 

qT 

qT + k L 

j (27r)2 Cp 

_ q\ + M'^ 

(kx+qT)2 + M^J 


(5.10) 

2 


where we have relabeled h = pq_i_ on the r.h.s.^® and kept the subleading term \n x ^ g, 
which is necessary to match to the diffusion equation. Indeed, one can check that, upon 
plugging Eq. (5.10) in Eqs. (5.2)-(5.3) and expanding consistently for x ^ g, the diffusion 
structure described in detail in App. D and in particular in Eq. (D.14) appears. The 
subtraction term then reads 


r^diff 

collsubtr 


[< 5/1 


T«/(p) + 



dSfjp) 1 (fdf{p) ' 

dp^ 2 d{p^Y 




coll 

subtr. 


(5.11) 


with 


dqi 


coll 


subtr. 


, 2f^r(p,a;i coll 

dujuj - - - 

dev 


Mu; 


Q^CrT ^ 


d^h 


g^CnCAT 


duj 


soft g Svrp^ 

d'^qx f d?k±_ Cpiqi) 


TT 




/ 


J (27r)2 

q± 


2h • ReF(h,p, w) 

k±+q± 


soft g 
2 


(2vr)2 7 (27r)2 Cp [q‘{+M^ (k_L+q±)2 + 


(5.12) 


where we have dropped the statistical factor on p. The subleading term in u in Eq. (5.10), 
while vanishing in Eq. (5.12), is critical in obtaining the necessary oc p~^ terms in 
Eq. (5.11). < T is a UV regulator for this region, as the approximations we have 

taken for the derivation of Eq. (5.12) fail when w ~ T. Indeed, there 5E becomes 
of the same size of f d^k±C(k±) and the LPM effect intervenes, so that the complete 
leading-order rate, as given by Eq. (5.3), is finite. 

This is a consequence of the form of Eqs. (2.3) and (5.1)-(5.2). In their derivation (see Eq. (2.6) in 
[21]) one integrates the effective 1 •<-/ 2 matrix elements over the transverse momenta of the final states, 
neglecting the small deviations from eikonality in the distribution functions, i.e. taking f{q') ~ fiq'p)- 
This makes the diffusion limit of Eqs. (5.1) and (5.2) insensitive to q, as it also happens when 5f is a 
function of p only, as we remarked at the end of Sec. 3.2. Transverse momentum broadening would enter 
in the diffusion limit of the collinear sector when taking the first correction to the eikonal approximation, 
which would take the form of q±V^f{q'p), assuming as usual p jj z, and would thus be suppressed by 
a factor of g^. Interestingly, this term would be responsible for the appearance of the double logarithm 
that has been recently pointed out in [17-19]. 

^®Due to the properties of the cross product, h and pqx have the same modulus but point in different 
directions, which is irrelevant in this case. 
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5.3.2 The conversion limit 


We now need to consider the q gq process with {p — u) ~ gT and the g qq one with 
either uj oi p — ui soft, which yields again a factor of 2. An altogether similar treatment 
then results in 




coll r] A 

subtr. UF 


coll 


subtr. 


(5.13) 


and similarly the antiquark and gluon terms have the same structure as their leading- 
order counterparts Eqs. (3.25)-(3.26), with the leading-order conversion rates replaced 
by subtraction rates (ih. These subtraction rates read 


r-nconv 


(p) 


coll 


subtr. 



. (5.14) 


.Qi + m'io (q± + k±)2 + m^_ 
Subleading corrections to the expansion of Eq. (5.8) are not needed in this case. 

5.3.3 The semi-collinear limit 

As we shall explain in more detail in 8, the semi-collinear regime refers to the region 
where q"^ ~ gT‘^ and no leg is soft, i.e. uj > T, p — uj > T. This in turn implies that 
h?/{uj{p — cj)) ~ gT"^ S> In this case, we have 

SE ^ SE{h,p,u) ^ (5.15) 

which is larger than f d^k±CF(k±) so that the integrated and symmetrized version of 
Eq. (5.8) becomes 


d?h 


2h • ReF(h) 


(fh f d?k± 2 k‘j_CFik±) 


semi—coll 


/ 


(27r)2 J (27r)2 5E2 Cp 




+^{p‘^ + (p-^f) 


. (5.16) 


Its equivalent for the g ^ qq process can be easily obtained. Since, as we shall show, 
the collision operator for the semi-collinear sector is conveniently formulated in the same 
form as Eq. (5.2), it suffices here to derive the subtraction rates, which read 


dT{p,Uj) collsubtr. 
duj I semi—coll 


g^Cn 

2'KP 


(1 ± n(a;))(l ± n{p — uj 


)) f — f 

”J 


X 
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d^kp k\CF{k±) 

(2vr)2 

Cpqi 

q^qg 


g^ qq 

>, 

g^ gg, 



(5.17) 
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where we have relabeled h? = p^q]_. 


6 The diffusion sector at NLO 


We now compute the NLO corrections to the diffusion coefficients of Eq. (3.23), 

Qnlo = Eq- (3.15) + 6q , Ql.nlo = Eq. (3.21) + . (6.1) 


The NLO corrections to q have been previously calculated, [33]: 


6q 


g^CnCAT^ 
327r2 T 


(37r2 + 10-41n2) , 


( 6 . 2 ) 


so we focus on the corrections to q^. These corrections will be the sum of three terms: 




SqL 


loop 


SqL 


coll 

subtr. 


SqL 


diff 

subtr. 


(6.3) 


where the three dq^ encode respectively the 0{g) loop corrections to longitudinal momen¬ 
tum diffusion, the collinear counterterm obtained in Eqs. (5.11)-(5.12) and a counterterm 
for a mistreated region in the LO calculation of q^. 


We start with dqi 


loop 


, the NLO soft contribution to Eq. (3.16) arising from adding 


one extra soft gluon. A first reduction in the number of relevant diagrams comes from 
the fact that, as observed in [29], we can write F~^ as F~^ = A~ — [H“,A+] and 

use the equation of motion of the Wilson line, D~j_U{x'^;Q) = 0, so that 


U{a- x+)[L>", A+(x+)][/(x+; h) 


(I7(a;x+) A+(a:+)C/(a:+;6)) , 


(6.4) 


i.e. the commutator acts as a total derivative {d~) and can be discarded in the dx'^ 
integration, provided that the boundary term vanishes. This is true in all non-singular 
gauges, where the A+ field vanishes at large x'^, such as the Coulomb or covariant gauge 
(or, trivially, in the singular = 0 gauge). Using translation invariance and shifting 
the integration by —x"*" the same trick can be applied to the other field strength insertion, 
so that in the end in Coulomb or covariant gauge we need to worry only about 


9 


r*+oo 


= — / dx^ Tt (U{—oo]x'^)d^ A {x'^ A (0)C/(0; —oo)), (6.5) 

J—oo 


where we have suppressed the trivial dependence of gauge fields and Wilson lines on 
the constant x~ and x± coordinates. A second simplification comes from noting that, 
similarly to leading order, at NLO operator ordering is not relevant in the soft sector 
in this case. Since in a first approximation ~ ~ Grr ~ Gi? ~ l/gGn, all gauge 
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fields must connect to the Wilson lines as r fields [29], so that we can replace the more 
complicated contour in Eq. (6.5) with a simpler adjoint Wilson line, i.e. 


Ql = 


9^Cr 


/ + 00 . 

dx+ ( 5 +^- “(x+)C/“^(x+; 0)9+A- ^(0) 

-00 ' 


J — OQ 

Its evaluation requires the computation of the diagrams shown in Fig. 8 .^® 


( 6 . 6 ) 



Figure 8 . Diagrams contributing to Sq^ 


loop 


at NLO. 


6qi 


diff 


subtr. 


arises from an 0{g) error we have committed in the previous determination 


of Ql to LO. Namely, we have used and resummed HTL self-energies in the LO calcu¬ 
lation, for instance, in Eq. (3.18), without worrying about the fact that the HTL loop 
integration extends down to zero momentum, where the hard approximations used to 
simplify the calculation of the HTL break down. In other words, the last two diagrams 
in Fig. 8 have already been included in our LO calculation, but using approximations 
which are invalid for small loop momentum. To fix this, we should subtract off the 
large-momentum limiting behavior of these diagrams when we evaluate them in the 
NLO computation. 

We present the details of the calculation of both terms in App. F. Here we just men¬ 
tion that the general structure corresponds to what was found for the soft contribution 
to the photon rate [27]. Schematically, the same sum-rule technology can be applied: 
the Wilson line propagators depend only on the minus components of the momenta, 
so that we can again deform the contour when integrating the plus component, which 
we call q^. This corresponds to expanding those diagrams for large, complex g+. The 
leading contribution should be of order ((?+)° and the subleading one of order ((?+)“^. 
Higher-order terms are suppressed and can be neglected. The leading, 0((g+)^) term, 
once integrated along the contour, will give rise to a linear divergence. An analogous 

coll 


linear divergence appears in dq^ 
divergences cancel. 


subtr. 


, as shown in Eq. (5.12). As expected, these linear 


^®3-point and 4-point vertices in these diagrams should be understood as including HTL corrections. 
However, after we deform the q~^ contour to large (complex) values, the contribution of the HTL vertices 
becomes small and they do not contribute to our final calculation. 
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A term behaving as 0{l/q^) at large has an interpretation of an asymptotic mass, 
which is why our LO result, Eq. (3.21), can be written in terms of the LO asymptotic 
mass. Therefore it is not surprising that the NLO correction is found by substituting 
the NLO form of the asymptotic mass + 6M^, as dehned in Eq. (5.6), into 

Eq. (3.21), and then expanding to linear order in SM^: 


Ml 

ql + Ml 


Ml + 5Ml _ Ml , ql 

ql + Ml + 5Ml - ql + Ml^^ {ql +Mlf 


(6.7) 


From Eq. (6.7) we thus obtain 




loop 


5qL 


coll 

subtr. 




diff 

subtr. 


T f qI^mI 

^ ^ J (2vr)2(gi+M^)2 

g^CnTSMl 1 ( 

4vr Ml J 



( 6 . 8 ) 


This simpleminded argument indeed reproduces the detailed explicit calculation of Ap¬ 
pendix F. 

Eq. (6.8) depends on a regulator As we will show, the dependence of this 

term on the regulator and the dependence of the semi-collinear region will cancel. This 
completes the evaluation of the diffusion sector to NLO. 


7 Conversion processes at NLO 


According to Eq. (4.3), the NLO corrections to the conversion sector take the following 


form: 



(7.1) 


where the sum is understood to give rise to the structure of Eqs. (3.24)-(3.26). The NLO 
conversion rates are composed of three parts, namely 


5r“4^(p) = drizi(p) 


loop 


^rni(p) 


coll 

subtr. 


^KTb(p) 


conv 

subtr. 


(7.2) 


The Hrst term on the right-hand side comes from the soft-gluon loop correction to the 
rates, as dehned by the Wilson-line operators (3.28)-(3.30). The second is the subtraction 
term from the collinear region, as obtained in Eq. (5.14), and the third subtracts the 
Hard Thermal Loop approximated leading-order calculation, result (3.31), in complete 
analogy to eiicountered in the previous section. 

The hrst and the third term can then be evaluated from Eqs. (3.28)-(3.30), by adding 
one extra soft gluon to the LO term in Eq. (3.31). A key observation is that the power¬ 
counting arguments that lead to the simplihed form for given by Eq. (6.6) apply 
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here as well: all soft gluons must connect to the Wilson line as r fields, so that their 
ordering is not relevant [29]. This implies that the fundamental and adjoint Wilson lines 
appearing in Eqs. (3.28)-(3.30) can be simplified to NLO to a simpler antifundamental 
line connecting the soft fermionic fields, i.e. 

p+oo 

(ix+(Tr[V’(x+,0,0_L)jit/F(0,0,0;x+,0,0_L)V’(0)]> , (7.3) 

r5“;(p) = ^r“"{p). (7.4) 

This corresponds to an effective abelianization of these operators, which, we note, are 
the same as those appearing in the fermionic sector of the Hard Thermal Loop action 
[31]. Indeed, an altogether similar abelianization happens for instance when obtaining 
the effective qqg HTL vertex in QCD. 

The diagrams necessary for the NLO evaluation of Eqs. (7.3) and (7.4) are shown in 
Eig. 9. However, we note that, in their abelianized forms, Eqs. (7.3) and (7.4) correspond. 





Figure 9. Diagrams obtained from Eq. (7.3) at NLO. The double line is the fundamental Wilson 
line and the black squares are the insertion of the soft fermion fields. 

up to the prefactors, to the soft-sector contribution to the NLO photon rate in [27]. 
Therefore we can directly use that result, which was obtained using the same sum- 
rule techniques employed in the previous section. Indeed, as we observed there, the two 
results are remarkably similar, the only difference being given by the different asymptotic 
masses. Here the relevant one is the quark one and we then have 

rpconv/ 4 ^ 9‘^Cf f (fq± q\ Sml^ ^ g'^CrSm 

P 4p J |'27r)2 Ibvrp 

where is given by Eq. (5.6) and the linear divergence in the collinear counterterm 
canceled an opposite one coming from the loop corrections. The logarithmic UV diver¬ 
gence has been treated with the same UV regulator used in the previous section 

for 6qL. 



(7.5) 


8 The semi-collinear region 

As we anticipated in Sec. 4, semi-collinear processes can be seen as 1 ■H' 2 splitting 
processes where the virtuality and correspondingly the opening angle are larger. Two 
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examples are drawn in Fig. 10. The scalings of this region are as follows: K ~ gT is soft 



Figure 10. Diagrams for two typical semi-collinear processes. In the first case the soft gluon 
is in the space-like Landau cut, whereas in the second case it is on its time-like plasmon pole, 
represented by the black blob. 

whereas the two hnal-state particles are quasi-collinear, i.e. with an increased virtuality 
and opening angle with respect to the collinear sector. The leading contribution then 
comes from T,q ~ gT, ql ~ gT"^, ^ gT ‘2 or, in the case of a democratic 

splitting, ~ E, q ^ gT, q‘j_ ~ gTE, Q'^ gTE. Naive power-counting arguments 
would suggest that the semi-collinear region should contribute to leading-order, as it 
is the largest slice of phase space where a soft gluon can attach to a 1 -H- 2 process. 
However, once all diagrams are summed and squared, a cancellation, discussed in [62] in 
the context of photon radiation, introduces an extra 0{g) suppression^*^. Furthermore, 
as we shall show, the contribution from time-like soft gluons, e.g. plasmons, is now 
allowed. 

The contribution ^( 7 ^ 61111-0011 collision operator can be written in the same way 

as the collinear one, as given by Eqs. (5.1)-(5.2), with the replacement of the collinear 
rates with semi-collinear ones. For instance, for quarks and antiquarks it reads 


( 8 . 


The derivation of the semi-collinear rates then requires the evaluation of processes of 
the form of Fig. 10, with p,q'^ S> q± ^ k±,k'^. Actually we have already evaluated 

^°The cancellation occurs because the transverse momentum of the split particles is larger than 
the disturbance from the scattering q\. In the limit that the disturbance is arbitrarily small, we would 
not expect it to induce a splitting. This cancellation can be seen at work in our derivation of Eq. (5.16) 
from Eq. (5.8). Upon enforcing semi-collinear kinematics on the latter, i.e. h pk^,a;kx, (p — iv)kx, 
all terms in square brackets vanish at first order in that expansion and only the next one gives a nonzero 
contribution. Without that cancellation the semi-collinear rate would indeed be leading-order. 
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these diagrams using the collinear expansion, since it is precisely these diagrams which 
give rise to the linear-in-collisions expressions we found in Subsec. 5.3.3. In particular, 
the subtraction term from the collinear region, Eq. (5.17), was derived by making an 
expansion in S> q±, and it still applies, under one condition. In evaluating the 
collision sector, we treated q_L ~ /c_i_ ~ gT, leading to 6E ~ g'^T. This let us neglect 
5E when working out the kinematics of the soft gluons, so that CR{k±) (see for instance 
Eq. (3.14)) is defined for k~ = 0 and hence only space-like gluons contribute to it. But if 
q'j_ ~ gT'^, q^ ~ T, or q‘j_ ~ gTE, ~ E, then 6E ~ gT and can no longer be neglected, 
opening up the time-like-gluon sector. In particular, when we put {P — Q) and {Q + K) 
on shell, we hnd, exactly as in the photon case [27], which we refer to for more details, 
that 


m + Kf) 


6{k- - SE) 

2^^ 


+ 0{^), 


5E 


PQ± 

2q+{p-q+) 


PQ± 

2u}{p — oj) ’ 


( 8 . 2 ) 


where the 0{^) correction comes from • k_L and always vanishes in the angular 
integrations. We see that 6E is exactly what we have used in Sec. 5.3.3. Therefore we 
must re-derive Eq. (5.17) with these somewhat different kinematics. A straightforward 
computation^^ shows that the findings in the case of photon radiation [27] generalize to 
the present case. Namely, the quantity 

/1^ . (8.3) 


physically interpreted as the transverse momentum diffusion coefficient and present in 
Eq. (5.17), should be replaced with its finite dE generalization. 


qm _ 

9^Cr 


d^K 

(27r)3 


6{k- - 6E) 


^lG--{K) + 2 G^t^{K) 


SE 




(8.4) 


which goes into Eq. (8.3) for SE —)• 0 and corresponds to the leading-order soft term (in 
Coulomb gauge) in the evaluation of the operator 


/*cio 

q^SE) = ^ (u^F/“(x+,0,0x)C/l'(x+,0,0x;0,0,0x)u^E^',(0)), 

J —oo 

(8.5) 

which was hrst introduced in the photon case [27]. In principle in the present case a 
more complicated “three-pole” operator should be needed [33]. However, at leading and 

^^Interestingly, the computation can also be performed using standard Soft Collinear Effective Theory 
(SCET) [63-68]. Indeed, we have Q = {q^,q~,q±) ~ A(l, A^,A) and K ~ A(A^,A^, A^), where A is the 
large scale, E or T, A <C 1 is the expansion parameter, either A ~ \/gTjE or A ~ y/g. These are then 
the standard scalings of SCETj. However, due to the cancellations mentioned in this Section, the 0(X) 
soft-collinear couplings [68-71] are necessary. 
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next-to-leading order it would reduce to a set of three two-body exchanges of the form 
of (8.5), with the appropriate Casimir factors [33]. 

Eq. (8.4) can be evaluated using Euclidean techniques, yielding [27]^^ 

qjSE) _ f 

9^Cr J (27r)2 


777-7 




26E^ 


(A:2 +6E^)(k‘^. +6E‘^ + ml) k^. + 6E^ 


( 8 . 6 ) 


However, since this momentum region has overlap with both the collinear and the hard 
regions, there are two subtractions which must be conducted, corresponding to the 
treatments already included in those leading-order calculations. Therefore we must 
compute the behavior of this momentum region under each of those limiting kinematics 
and subtract them. The collinear case is treated by subtracting q from q(6E). For 
the hard region, we take the soft, bare limit (?7s(/c°) —>• T/k^, p{K) —)• p^^\K) = 
2tt sgTi{k'^)5{K'^)) of Eq. (8.4), yielding 
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The full semi-collinear rate is then obtained by replacing 


q 

9^Cr 


qjSE) 

9^Cr 


q 

9^Cr 


q{m 

9^Cr 


hard 
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in Eq. (5.17), which yields 


dT{p,uj) 


du 


9^CrT, 


d‘^kj 1 
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X 




^E^Ca[ 1 - X + x"^] 9^99 
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7777 


{kl + 6E‘^){kl + 6E‘^+ml) kl + ml 


( 8 . 8 ) 


We stress that the collision operator has the same form as Eqs. (5.1)-(5.2). 

The q± integration in Eq. (8.8) is to be understood as IR-regulated by pJ^^. In 
App. G we show how the small-w-and-gx region gives rise to IR logarithms that cancel 
the dependence of the diffusion and conversion sectors. We also give some details 

of how the transverse integrations can be carried out analytically. The w integration 
remains to be performed numerically. 

^^The Euclidean evaluation combines the time-like (plasmon) and space-like (scattering) contributions. 
From Eq. (8.4) it follows that, once k~ is integrated over the li-function, plasmons contribute for > 
k\/{25E), while space-like gluons contribute for k'^ < k\/(25E). In order to disentangle the two 
contributions one would have to proceed numerically. 
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9 Summary and conclusions 


The main aim of this paper has been to show how the propagation of highly energetic 
quarks and gluons through the QGP can be described at leading- and next-to-leading 
order by a Boltzmann equation encoding the interaction between these hard particles 
and the thermal and soft constituents of the plasma. Sec. 2 has been devoted to a brief 
review of the LO kinetic approach introduced in [21] and implemented in MARTINI. 
As Eq. (2.1) summarizes, the two processes it incorporates are 2 •(-)• 2 scatterings with 
the thermal medium constituents and 1 •(-)• 2 collinear splittings induced by the soft 
background. 

In Sec. 3 we have shown how this approach is not optimal beyond leading order, 
where the distinction between the two classes would blur and the resummed matrix- 
element approach to 2 2 scattering would become cumbersome. With these motiva¬ 

tions, we have reorganized the LO collision operator into four separate processes which 
provide a sufficient description at NLO. They are large-angle scatterings, i.e. 2 -H- 2 scat¬ 
terings with 0(1) angles or equivalently 0{T) or larger transferred momentum, diffusion 
processes, caused by soft gluon exchanges, which preserve the identity of the hard par¬ 
ticles while slightly affecting their momentum, conversion processes which instead turn 
quarks into gluons and vice versa through soft quark exchange and finally collinear 
processes, corresponding at LO to 1 •(->• 2 processes. In Sec. 3.1 we described in detail 
our description of large-angle processes, which require regularization to be kept sepa¬ 
rate from diffusion and conversion ones, as shown in Eqs. (3.5) and (3.7). Sec. 3.2 has 
been dedicated to diffusion processes, which are described by an effective Fokker-Planck 
equation, Eq. (3.8). The three physical effects of drag (energy loss), longitudinal and 
transverse momentum broadening are encoded in three corresponding coefficients. The 
requirements that the Fokker-Planck picture be equivalent to the Boltzmann one and 
that it approach equilibrium can be used to write the drag coefficient in terms of the other 
two, as per Eq. (3.22). The two momentum diffusion coefficients can then effectively be 
described by field strength correlators along Wilson lines on the light-cone direction of 
propagation of the hard particle, as in Eq. (3.16). The calculation of the transverse 
momentum diffusion coefficient q is mapped to a Euclidean one [33], whereas for the 
longitudinal momentum diffusion coefficient we introduce a sum rule which, through 
the analytical properties of amplitudes at light-like separations, makes it sensitive only 
to the gluon dispersion relation close to the light cone (see Eq. (3.21)). Similarly, con¬ 
version processes are shown in Sec. 3.3 to be described by effective Wilson line operators 
(Eqs. (3.28) and (3.29)), which are also computed through an equivalent light-cone sum 
rule mapping them to the quark dispersion relation. The UV log-divergence of the diffu¬ 
sion and conversion processes cancels with the opposite IR one in large-angle scatterings. 

In Sec. 4 we introduced the NLO extension of this reorganized approach. All pro¬ 
cesses, with the exception of large-angle scatterings, are sensitive to 0{g) corrections 
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arising from the interactions with the soft background. Furthermore, some care is nec¬ 
essary in avoiding double countings in slices of the phase space, which were included at 
LO, where some particles become soft, introducing the need for a set of subtractions. 
The remaining sections are then devoted to the details of each process at NLO. In Sec. 5 
we discuss the collinear region, which is sensitive to 0{g) corrections in the interactions 
with the soft background that induce the splitting, as well as in the dispersion relation 
of the hard and thermal particles. We further identify all necessary subtractions. 

Sec. 6 is dedicated to 0{g) corrections to diffusion. In treating q, we employ the NLO 
determination of [33], whereas for we perform the calculation using the light-cone sum 
rules introduced before. The details are to be found in App. F. The result is surprisingly 
simple: it just amounts to considering the soft correction to the gluon dispersion relation 
close to the light-cone (see Eq. (6.7)). Similarly, conversion processes are dealt with 
using the fermionic analogue of the same sum rule and require the inclusion of the soft 
correction to the quark asymptotic mass, as in Eq. (7.5). Both and the conversion 
rate at NLO show an UV logarithmic divergence, which is removed once a new process, 
which only starts to contribute at NLO, is considered, the semi-collinear process. As 
illustrated in Sec. 8, this process appears as a bridge between the diffusion/conversion 
sector on one side and the collinear on the other. Indeed, while retaining a collinear 
kinematics, it shows relaxed constraints, going beyond strict collinearity and allowing 
the interactions with the soft background to be not just space-like (soft scatterings) but 
also time-like (plasmon absorption/emission). Eor its evaluation a modified form of q, 
q{6E), is introduced in Eq. (8.5). It accounts for the changes in the small light-cone 
component p~, which are no longer negligible. Euclidean techniques are used for its 
computation, as per Eq. (8.6). 

We would like to emphasize the importance of Euclidean techniques, which map the 
calculation of C(A:_i_), q, q{SE), SM^ and into simpler calculations in dimensionally- 
reduced EQCD. Similarly, light-cone sum rules reduce the computation of q^ and of 
the conversion rates to the determination of the gluon and quark asymptotic masses 
at leading- and next-to-leading order. Without these recent theoretical developments, 
rooted in the causal properties of amplitudes at light-like separations, the calculations 
presented here would have required extensive, cumbersome numerical integrations over 
the intricate structures of loops composed of HTL propagators and vertices. Eurther- 
more, as we have mentioned. Euclidean techniques also allow lattice determinations. 
The first measurements of q and C(x_l) have recently been reported [51, 52], opening up 
a new avenue of research. All other Euclidean operators can be computed on the lattice 
in the same way, creating the tantalizing possibility of a factorized approach to kinetics, 
where perturbation theory is used at the thermal and hard scales to compute the large- 
angle scatterings and the splittings, whereas the 3D lattice is employed at the soft (and 
ultrasoft) scale to determine non-perturbatively the transverse diffusion processes and 
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the scatterings leading to collinear radiation. 

A very important point we have not addressed in this paper, leaving it to future 
work, is the impact of the NLO corrections we have introduced on calculations of jet 
modification and their comparison to experimental data. As we mentioned, the Monte 
Carlo event generator MARTINI implements a kinetic approach corresponding to the one 
described in Sec. 2. This makes it an ideal candidate for the inclusion of the NLO cor¬ 
rections. Indeed, the reorganization of the LO collision operator in terms of large-angle, 
diffusion, conversion and collinear processes is underway, as well as the implementation 
of the NLO corrections. This could also be easily complemented by the inclusion of non- 
perturbative input, such as the existing determination of q and, should they become 
available, future determinations of q{6E) and of the asymptotic masses. It would also 
be interesting to study the angular structures of jets with this numerical implementation 
and compare it with the recent order-of-magnitude perturbative estimates from [72]. 

We remark that it is difficult for us to gauge a priori the impact of NLO corrections 
relative to LO. The recent NLO calculations of the thermal photon [27] and low-mass 
dilepton [28] rates, which include many of the features presented here, such as Euclidean 
techniques, light-cone sum rules, semi-collinear and collinear processes, showed how the 
NLO corrections naturally grouped into two classes of large, and largely canceling, con¬ 
tributions. The positive corrections were due to NLO modifications to the collinear pro¬ 
cesses, caused by the increased soft scattering rate and the reduced asymptotic masses, 
while semi-collinear and conversion processes decreased the rate by a similar magnitude. 
The large cancellation between these contribution is mostly accidental and furthermore 
depends significantly on the details of the medium, such as the numbers of colors and 
flavors. So, while we anticipate similar cancellations for the present energy loss case, we 
are at present unable to quantify their impact in more detail. 

Finally, we believe that the approach presented here should go much of the way 
towards making possible NLO kinetic theory calculations of the shear viscosity and other 
transport coefficients of QCD. However, we have not resolved the issue of keeping track of 
where the energy from a soft scattering shows up amongst the other (thermal) particles, 
which so far prevents us from a true NLO calculation of QCD transport coefficients. We 
hope to return to this issue in the future. We do note, however, that for cases where the 
momentum dependence of the off-equilibrium distributions is isotropic, such as studies 
of isotropic thermalization [73, 74], an extension to NLO appears within reach. 
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A Notation 

We now summarize our notation. We will use capital letters for four-vectors, lowercase 
italic letters for the modulus of the spatial three-vectors, and the mostly-plus metric 
= Diag [—h++], so that - pg- 

For convenience we will mostly work in the Keldysh, or r, o, basis of the real-time 
formalism for the computation of thermal expectation values. The two elements of 
this basis are defined as cjir = {4>i -|- 02 )/ 2 , (f^a = 4>i ~ 4>2, 4^ being a generic field and 
the subscripts 1 and 2 labeling the time-ordered and anti-time-ordered branches of the 
Schwinger-Keldysh contour respectively. The propagator is a 2 x 2 matrix, where one 
entry is always zero and only one entry depends on the thermal distribution, i.e., 

^ f Drr Dra ^ ^ n(p°)) {Dr - Da) Dr\ ^ 

\Dar Daa) Da 0 ^ ^ ^ 

where Dr and Da are the retarded and advanced propagators, the plus (minus) sign 
refers to bosons (fermions). n{p^) is the corresponding thermal distribution, either 
= (exp(p*^/T) — 1)“^ for bosons or np{p^) = {exp{p^/T) + 1)“^ for fermions. 
We also define the spectral function as the difference of the retarded and advanced 
propagators, p = Dr — Da- We will denote the gluon propagator by G and the quark 
one S. 

We will adopt strict Coulomb gauge throughout. The treatment of soft momenta 
in propagators and vertices requires the use of Hard Thermal Loop (HTL) resummation 
[30]. For convenience we list the Coulomb gauge retarded HTL resummed propagators 
for fermions and gluons in the next section. 


B Hard Thermal Loop propagators 


In this section we detail our conventions for the HTL propagators. Fermion propagators 
are most easily written in terms of components with positive and negative chirality-to- 
helicity ratio. The retarded fermion propagator reads 

SR{p) = hpUP) + KSRiP). (B.i) 

where 


S^P) = 


pO zp (p-p S±(p0/p)) 


p°=F 


P + 


m: 


2p 


1-^ln 

2p 


p° +p 

p^ — p 


p0_p0_|_2e 

(B.2) 
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where the upper (lower) sign refers to the positive (negative) chirality-to-helicity compo¬ 
nent. The projectors are hp = ■ p)/2. Here is the fermionic asymptotic mass 

squared, defined such that the large-momentum dispersion relation for helicity=chirality 
fermions is Pq = p'^ + We similarly define the asymptotic gluonic mass M^. At 
leading order, their values are 


= — 2 . 
oo 2 


^2 

6 


iVe + 


N 


2n-<2 




= 2m^ = Cp 


g^T 


(B.3) 


where we have also shown the relations to the more commonly used Debye mass rrio and 
quark “mass” niq. 

Gluons are described in the strict Coulomb gauge by 


GS'(Q) = 


2 I 2 

q + nip 


1 (f ^ + q + te 

I-In - 

2q q^ — q + ie 


GliQ) = ( 5 *^ - <fq=)Gl{Q) = 


_ qiqi) 


to-f- Vi (|| 


q^ J 2q q^—q 


(B.4) 


qO=qOpie 

(B.5) 


The other components of the propagators in the r, a basis can be obtained through 
Eq. (A.l). 


C Longitudinal momentum diffusion from Wilson lines 

Eq. (3.16) is based on eikonalization, which naturally happens since p is considered in¬ 
finitely larger than all other scales at leading order. As such, it can be easily verified 
that the pertnrbative expansion of Eq. (3.16) agrees with the rate-based definition (3.11) 
at leading and next-to-leading order. We believe that, in the presence of a consistent 
UV regulator^^ Eq. (3.16) is correct to all orders in g at the leading order in l/p, up 
to possible Wilson lines along the x~ direction at x"*" = —oo. Indeed, we believe that 
Eq. (3.16) can be rigorously obtained in dimensional regularization using SCET, analo¬ 
gously to what has been done in [44, 75] for q. We sketch here a simplistic derivation. 
Since we are interested in the differential-in-q^ rate for a fast particle propagating with 
pO _ j|- jg natural to expect from the eikonal approximation a correlator of the form 

(27r)^ = (27r)^=^lin^l J ^(Tr[/(-L/2, L/2; x-)C/(L/2,-L/2; 0)>, 

use a UV cutoff fj,± in this paper; for a more rigorous treatment we could use dimensional 
regularization, or the introduction of a mass with the limit p —> oo taken holding m/p small but finite, 
which produces a “dead cone” which renders radiative effects finite. 
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where we have used the fact that in the infinite-p limit and for simplicity we 

have introduced 


[/(a"*", 6’*“; c ) = Pexp 



and for further convenience 


U{a- 


b ; c^) = P exp 




(C.2) 


(C.3) 


The Wilson line at = 0 is supported on the time-ordered branch of the Schwinger- 
Keldysh contour and conversely the other one is supported on the anti-time ordered 
branch, corresponding to the amplitude and conjugate amplitude entering the definition 
of the rate. Indeed, so far the techniques used in [44, 75] are exactly applicable here as 
well, so that Eq. (C.l) is also formally justified within SCET. It is however not gauge- 
invariant. Eollowing the steps of [44], we conjecture this form for its gauge-invariant 
dressing; 



(L/2,x-) 


(i/2,0) 


(L/2,-oo) 


Figure 11. The Wilson loop giving rise to qi,. The horizontal axis is the -I- axis and the vertical 
one is the — one. The points are given in coordinates, the constant transverse one is 

not shown. 

(27r)—— = lim — [ —(TrI7(—oo, —T/2)[/(—L/2, L/2; x“) 

dq+ L^oc L J dji 

xU{x-, -oo; L/2)U{-oo, 0; L/2)[/(L/2, -L/2; 0)f7(0, -oo; -L/2)). (C.4) 

The operator defined by Eq. (C.4) is sketched in Fig. 11. This particular ordering 
corresponds to having the upper three connected Wilson lines on the anti-time ordered 
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branch of the Schwinger-Keldysh contour and the lower three on the time-ordered one. 
The “handle” on the bottom right corner can be trivially annihilated, but the same is 
not true for the one at the bottom left, since time-like separated fields appear between 
the two vertical Wilson lines there. 

Finally, by using the dehnition of and convoluting Eq. (C.4) with the latter 

can be replaced by derivatives which, when acting on the Wilson loop, introduce the 
F~^ electric fields. Once the q'^ integration is taken (with infinite cutoff, hence the strict 
validity in dimensional regularization only), the Wilson line operator (C.4) is squeezed 
to the form of Eq. (3.16) plus a surviving “handle” along x~ at x"*" = —oo. This handle 
is irrelevant in non-singular gauges and even in the light-cone gauge A~ = 0 it can be 
neglected at LO and NLO.^^ The same would not be true for dp^/dt (in the p —)• oo 
limit), where we would encounter a single F~^ insertion (at x'^) and the handle would 
be critical in obtaining a gauge-invariant leading-order result. 


D Leading-order matching 


In this section we shall prove how the diffusion-|-conversion-|-large-angle scattering is 
equivalent to the dressed 2 2 processes of [21]. 


D.l Diffusion matching 


For simplicity, we only consider diffusion matching for 
previously illustrated in Eq. (3.2). In the prescription of 
using the identity 


s^+u^ 1 1 (s—u)^ 

^ 2^2 F ’ 


the qiq 2 O- qiq 2 contribution 
[21], that process is treated by 

(D.l) 


and the replacement 


G%{P-P') {P+P'Y {K+K') 


(D.2) 


where G^{Q) is the retarded HTL propagator, as given by Eqs. (B.4) and (B.5). Upon 
plugging this into Eq. (3.5), putting the IR regulator to zero, introducing instead an UV 

In the A~ = 0 gauge the leading-order term arises from the < 9“T+(a;+)9“T'''(0) > propagator. At 
NLO only soft gluon corrections to that propagator can contribute. Soft gluons connecting the propagator 
and the handle cannot contribute: as we have remarked, in the soft limit ~ G"' ~ G^’^ ~ XjgGn, so 
that these soft gluons have to connect to the handle as r fields and their contribution cancels between 
the two branches of the handle. 
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regulator T S> S> gT and consistently expanding for oj, q± rsj gT, we get 


/*+oo /*^9_L n 

dioj^ dknAk)[l-n^{k)]x (D.3) 

k' ^ 


2\G^n{Q)f+^^\Gl{Q)fy 

loT — co‘^{l—2np{k)) 


k + u} { 1 — 


P 


+ 0{cv^ql)\x 


2^2 


-5f^{p) + 


uT - u^{l-npik)) d6P^ip) d6f^^{p) 


dp 


dp^ 


where G{Q) is understood to be G{LJ,q = yu;"^ + cp^) and we have omitted the stimu¬ 
lation factor and its derivatives on the outgoing hard leg, as they are all exponentially 
suppressed. Up to higher-order corrections we can put the lower integration limit for the 
k integration to zero^®. The terms within the first set of curly brackets come from the 
expansion of the (/>-averaged matrix element, whereas those in the second set come from 
the expansion of the distribution functions. We remark that the square moduli of the 
propagators on the first line are even functions of oj. Hence, the terms that would naively 
be of leading order in this expansion in g, i.e. those multiplying 6f and its hrst derivative 
on the third line, vanish in the integration. Keeping only the surviving, even-in-w pieces 
and performing the k integration we have 


D - 


g^T^ Cpg^ 

6 327r2 


P + OO 

/ dio 

J—oo J 0 


dq±q± 


p-{Q) 




p dp 


1 + ^) 

P J 


+ T 


dSf^^jp) 

dp^ 


(D.4) 


where we have also used the following relations, based on the explicit form of the prop¬ 
agator in Eqs. (B.4) and (B.5): 

/(Q) = G^(Q) - Gi(g) = irmly |G^|' , 

p^{Q) = GliQ) - GliQ) = irml^ (l - \Glf . (D.5) 

We also remark that the next terms in the soft expansion, i.e. those 0{uj'^,qj), which 
naively would contribute to relative 0{g), give rise again to a vanishing odd integration 
and thus contribute only to 0{g‘^). For this reason they can be neglected.^® Similarly, the 

The same is not possible when there are bosonic degrees of freedom associated with k, due to Bose 
enhancements. There, one needs to consider this region with care; this region is part of the semi-collinear 
processes. 

^®In order to obtain the explicit form of these 0{LtT ,q\) terms, the prescription illustrated in Eqs. (D.l) 
and Eq. (D.2) is no longer sufficient, as Hard Thermal Loops need to be included also on less IR-sensitive 
terms. However, on general grounds, the expansion can only give rise to even powers of cj at that order, 
as we have checked explicitly. 
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next order in the expansion of the distribution function is also odd and vanishes. Hence, 
a genuine 0(g) correction can only arise from adding soft gluons to these diagrams. 

An analogous expression can be obtained in the case of a p-dependent 6f. It reads 




_ p 

dq± q± 




^ 2T\ r.2 

1H I q± 

P J P . 


(D. 6 ) 


Now let us look at the coefficients entering Eq. (3.8), as defined in Eqs.(3.10), (3.11) 
and (3.12). The differential rates appearing there can be easily inferred from the loss 
term of the collision operator (3.2). Applying the same steps that led to Eq. (D.4) we 
have that the contribution from scattering with a quark q 2 to dp^ / dt for a quark reads 


dpL 

dt 


91 



(2|G|j(Q)|V^|G^(Q)|') 


xk 

r, f 


k + uj [ 

1-- 


[ V 

pj \ 


+ 0 {u:\ql) 



{k) [1 - Hpik + w)]. 


(D.7) 


where we have for clarity left the stimulation factor [l — np{k + oj)\ unexpanded in a; ~ 5 . 
Using Eq. (3.6), i.e. q^ = u + q\/{2p), one sees again that the naive leading order in g 
leads to a vanishing a; integration, which is at the base of the Einstein relation relating 
dpp/dt and q^ in the p —)• 00 limit. The leading, even-in-cu terms then yield, upon 
performing the k integration and using again Eq. (D.5) 


dpL 

dt 


91 


Cpg^ 

6 327r2 



dq^q± 


P (Q) 

m%uj 



(D. 8 ) 


where we have again not considered the 0{g‘^) correction from the 0 (ai^, q^) terms in the 
expansion of the matrix elements and from the expansion of the distribution functions. 

Returning to q^^ it is immediate to see that, at leading order, only the term in 
(g^)^ contributes and other terms are actually suppressed by a factor of g^. Hence q^ 
reads at LO 


qL 
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91 


g^T^ Cpg^ 
6 lOvr^ 


r>+00 


doj 


Jo 


dq± q± oj 


2 P iQ) 

nn%uj 


(D.9) 


Similarly, for q one has q]_ = qj_, up to odd corrections or 0{g^) terms, so that one 
obtains the well known result (see Eq. (3.15)) 


g^T^ Cpg^ 
6 IGtt^ 


_ 2P-{Q) 

duj / dq±q±qj_ - 5 -, 

-00 Jo 


(D.IO) 
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We furthermore remark that in and q corrections in 1/p enter only at 0{g^), due again 
to the 0{g) term being odd in oj. For this same reason the evaluation in (w,gj_) and 
iQi) coordinates is equivalent. This justifies our evaluation of the (p-independent) 
NLO corrections to and q in the latter coordinate set. 

Let us obtain the complete leading-order dp^^/dt, q^ and q. To this end, one has 
— 1) quarks and antiquarks that are distinguishable from qi, and hence 4(Nf — 1) 
contributions in the form of Eqs. (D.8), (D.9) and (D.IO), the extra factor of 2 coming 
from the sum over hnal states (see footnote 5). The contribution from qiqi scattering 
accounts for two times those equations, as the u-channel contribution is identical, and the 
qiqi accounts for another two due again to hnal state symmetries. Hence the contribution 
from all quark scatterings account for a factor of 4:Nf. Using altogether similar steps 
one can show that the contribution from qig scatterings amounts to a factor of 8Nc, 
that the complete leading-order expressions are 


so 
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dt 
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(D.ll) 

(D.12) 

(D.13) 


which agree with Eqs. (3.21) and (3.15). In the case where the hard particle is a gluon, 
one obtains the same expressions with Cp replaced by Ca- 

We can now see explicitly, by comparing Eqs. (D.ll), (D.12) and (D.13), that 
the equilibration condition Eq. (3.22) is obeyed at leading order. Finally, let us take 
Eq. (3.8), and substitute the equilibration condition. This yields, for a p-dependent 6f, 

d^5f{p) 


C^^[5f] = -^6f{p) - 
a I Ji Tp-’^^’ dp 


qL 


1 1 
p ^ ^ 


1 




dp^ 


(D.ld) 


which matches with the structure of Eq. (D.d). In the p-dependent case we recover 
instead Eq. (3.23), which also matches with Eq. (D.6). We have thus explicitly shown 
how the effective diffusion picture of Eq. (3.8) matches exactly at leading order with the 
standard treatment of dressed matrix elements. 


D.2 Conversion matching 

Let us consider more in detail the t-channel quark exchange to the qiqi -H- gg process, 
as introduced in Eq. (3.7). The resummation of HTLs in the t propagator, as per the 
prescription of [21, 39], gives 

# u 

Jo 


2tt t 
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{w-qy Sj(g)Sj((3) + (!.. + qY SpQ)SAQ) 


i + o|“ 

T p 


(D.15) 
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where we have used our parameterization (B.l) of the quark propagator. Although not 
immediately obvious in a naive expansion, all 0{g) corrections do take the form of an 
odd function of w, whereas the leading order, i.e. the terms in square brackets, are even. 
By using the explicit form of the propagator in Eq. (B.2) the expression above simplifies 
to 


d(p u 

Jo 27r7 


2pk 


mt 


[(g - uj) p+{Q) + {q + u) p-{Q)] 


p 


(D.16) 


where p±{Q) = S^{Q) — S^{Q). Similarly the expansion of the statistical factors, as in 
Eq. (3.27), leads to odd terms in a; as the only possible 0{g) corrections. 

Hence, summing all contributions^^, the conversion part of the collision operator for 
a quark i reads 


c^rw] 


g^c] 


,4r^2 
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STT^m^p 


- 5/^(p)} j dkknp{k)[l + nB{k)] j 




p+(Q)(l--)+p_(Q)(l + - 




dq_i_q± 

(D.17) 


and corrections are naturally suppressed by because of the even ui integration. Car¬ 
rying out the k integration leads to 


c^^rw] 


Cf 
IGvr^ p 


(*00 P j-^q 

|5/''i(p) - 5/3(p)} J duj dq±q± 




p+(Q)(l--)+p_(Q)(l + - 


{^ + o{g^)), 


(D.18) 


Finally, the oj integration can be performed using the sum rule in [27, 58], leading to 

C“"'’l'5/l = f—Uf'“{p)-Sf‘'{p)} r" 2 ■ (D-“) 

^ Svr p I ) Jq q^^+ 

which matches with Eqs. (3.24) and (3.31). 

The conversion operator for gluons can be easily checked using the same approach. 
The case of a p-dependent distribution function is also a straightforward generalization. 
It too matches with the results of Sec. 3.3. 


E Solving the integral equation in position space at LO and NLO 

The most convenient way to solve Eq. (5.4)^® is by Fourier transforming h and q_L into 
impact-parameter variables, as first proposed in [45]. In this way the convolution over 

^^This amounts to the u/t and t/u terms for q\qi -f-j- gg, as well as the s/u one for qig -f-j- qig. The 
non-underlined u/s there can be easily shown not to contribute at leading and next-to-leading order in 
g, whereas the t/u and s/u terms become identical to Eq. (D.16). 

^®The g <rp>- qq case is again not dealt with explicitly. 
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the collision kernel C(/c_l) diagonalizes, turning an integral equation into a differential 
equation. Furthermore, the source on the left-hand side becomes a boundary condition 
at b = 0 and the desired final integral, Eq. (5.3), becomes a boundary value of the ODE 
solution. Specifically, defining 


we have 

/ i‘2, 'L^ 

^2h.F(h) = Im(2Vfe-F(0)), 

and Eq. (5.4) becomes 


(E.l) 


(E.2) 


-2fV(52(b) = 


2puj{p — Lo) 


(p(p - a;)m^^ -kpw - uj(p - uj)m^p - V^) F(b) 


+ (cMkI ‘-I'-) ) F(b). (E,3) 


with 

I |^(l-«'"'’‘‘")c«(U). (E.4) 

As we mentioned in Sec. 5.2, for generic kinematics 0{g) corrections enter then in 
two places: the effective thermal masses squared and the collision kernel C{k±) get 
0{g) corrections which modify Eq. (E.3), 

^oop,LO-l-NLO ^oop ^^oopi (^•^) 

l^/?LO-l-NLo(^) = ^Ri^) + (E.6) 


The NLO thermal masses have been given in Eq. (5.6). The NLO collision kernel is 
computed in [33] in momentum space; the Fourier transformation into impact parameter 
space has been performed in [27]. The expressions are sufficiently cumbersome that we 
have decided not to repeat them here. In [33] it was also explicitly shown that “three- 
pole” contributions are absent at NLO, so that the sum of two-body (dipole) interactions 
on the second line of Eq. (E.3) still holds. 

Eq. (E.3) is then solved perturbatively, by treating F(b) formally as an expansion 
in powers of Srrioo, bC\ F(b) = Fo(b) -|- Fi(b) and expanding to first order. The 

zero-order expression is just Eq. (E.3), while at the linear order the expression reads 
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(E.7) 
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where the leading order solution Fo(b) acts as a source term in the differential equation 
for Fi(b). We refer to [27, 28, 76] for details on the boundary conditions and the 
numerical evaluation of these equations. 


F Longitudinal momentum diffusion at NLO 


In this appendix we present the details of the calculation of to NLO. We will not 
explicitly consider diagrams with HTL vertices: as in the photon case, their contribution 
can be shown to be suppressed once the contour is deformed away from the real axis. 
Furthermore, when performing such deformations, we will not explicitly keep track of 
contributions from certain causality-violating poles at q~^ = q~ j2± iq± {(f' = 0), which 
are artifacts of our gauge choice and cancel in the final sum over diagrams, as they must 
[27]. 


F.l The rainbow diagram 

Let us first go through the diagrams shown in Fig. 8, contributing to dq^ 


loop 


(we will 


drop the loop label to avoid clutter). We will label hrst “rainbow” diagram, shown in 



Figure 12. The rainbow diagram 
Fig. 12, r. Its contribution reads 


5qi 


= -g^Cn 


/ +00 pa 

dx^ / 
-oo J 0 


r* + 00 PX~^ PX^' P ^4^ 

dx-l dx-'l dx-"j^^ 


—oo Jo 


d^K 


(F.l) 


where, as we remarked in Sec. 6, the specific ordering of the two propagators is not rele¬ 
vant to NLO, as long as they receive a Bose enhancement. The Wilson line integrations 
yield 


dqL 


r 


q^GrCa 



_ Kg""? _ 

{q~ + ie)‘^{q~ + k~ + ie) 


adv Gj.^ (Q)G'rr (-^)> 


(F.2) 
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where “adv” stands for the advanced —ie —)• +ie counterpart of the first term in round 
brackets. As in Sec. 3.2 and Fig. 7, we set out to perform the integration in the com¬ 
plex plane. The integral is very sensitive to large g"*" due to the in the numerator; 

but, contrary to the leading-order case, q~ is not fixed to be zero; also, G~~{Q) con¬ 
tains the statistical function ng{q^) ci^T/q^ = T/[q'^ +q~/2). Applying some numerator 
algebra to these terms, we obtain 


T{q^? ^ ^ , nq-f 

g+ + g-/2 2 4(g+ + g-/2)- 


(F.3) 


The hrst term yields the contour deformation, the second will vanish as we shall show 
(no poles and no contour contributions) and the third can be dealt with using Euclidean 
technology. 

We start with the contribution from the first term, with additional label (a) for 
arc. Upon deforming g"*" away from the real axis, the retarded propagator turns into 
Eq. (3.20), so that 


5qL 


(a) 


g^CRCAT [ 
Jc 

xG--{K)^ 


d^Q f d^K 


(27r)^ y (2vr)^ \(g +ieY{q +k + ie) 


1 + 


2g+g- - Ml 


q+J 2q+q- - q\- Ml 


Cr 


+ Ca: 



(E.4) 


where Cr and Ca are the retarded and advanced deformed contours, as defined in Sec. 3.2. 
The contribution from the latter is not shown explicitly. Let us define and, for later 
convenience, 5Eq_|_k and 5E'q_k as 


= 2,+ ^ 


^-®q+k = 


(q± + k^)^ + Ml 


2 

oo 


2g+ 


<5Eq„k = 


(q L - k L)^ + M, 


2 

oo 


2g+ 


(E.5) 


When deforming above the g"*" axis we can then close the q~ contour in the lower half¬ 
plane, picking the pinched q~ = (JEq retarded pole (and conversely for the Ca contribu¬ 
tion). This yields 


dqL 


(a) 


Mq^GrCaT 

xG--(A)4 

g+ 


r dq+d‘^q_L 
'cg (2vr)3 





+ Ca- 


dE^iSE^ + k- - ie) 


(E.6) 


The final expression, up to order l/g'*' terms, reads 


5qi 


(a) 


= q^GrGaT 


'Cr 


dq+ f d‘^qs_ f d^K qlG-E{K) 
(27r)4 2{q+)^El 


2-n 


f d^qi. f 

J (27r)2 J 


idEct 

Trd{k ) -j- _ ^ 


{k~ — ieY 


ECa , 
(E.7) 
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where we have used the symmetries of the integrand to express the leading-order term 
as a (5-function of k~. 

We now inspect the second term, labeled (s) 


(-) g^GnGA 

f d^Q 

fd‘K 

^ ' adv^ 

r 2 j 

(27r)4 j 

(27r)4 

V{q-+ie){q-+k-+ie) J 


Sql 

When deforming on Cr and Ca we have 

(«) g^CnCA f dq+ f dq-£q± f d^K 


(F.8) 


5qi 


L 


/ 


Ch (27r) 7 (27r)3 J {2tt)^ \{q - ie){q +k - ie) 


xG--{K)- 


1 + ^ 


- Ml 


{q+Y \ q+J 2q+{q-- 5E^ + ie) 
The q~ integration can be performed as before, yielding 


+ Ca- 


(F.9) 


dqi 




= -l- 


.g^CnCA f dq+ f d^q^ f d^K 


L 


/ 


(27r) J (27r)2 J {2ttY \ {5E^){6E^ + k-- ie) 


xG--(K) 






(F.IO) 


which goes like l/{q'^)‘^ and hence is irrelevant. This can be easily understood by noting 
that the pinched poles in q~ force q~ ~ !/<?''', so that the factor of q~/q~^ of this term 
with respect to Eq. (F.4) behaves like l/{q'^)'^. 

Finally, we look at the Euclidean term, labeled (e), which reads 


6qi 


(e) 


= g^CnCA 


d^Q 


^27r%- + k-)lG--{Q)G--{K). (E.ll) 


(27r)4 J (27r)4 M 

This term will be canceled by an opposite term in another diagram. 


F.2 The crossed rainbow diagram 



Figure 13. The crossed rainbow diagram 
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The amplitude of this diagram, shown in Fig. 13 and labeled + for cross, reads 


5qL 


= +g^CRCA 


/ +c» pa 

dx'^ / 
-OO J 0 


r + ' 


dx~^' f dx^" 

Jo 


I 


d^Q 


d^K 


{2ttY J {2ttY 




(F.12) 

The sign is opposite to Eq. (F.l) because of the different ordering of the color matrices. 
Doing the x~^ integrations we obtain 

d^Q r d^K 


6qr 


I 


+ 9 R A J J (27r)"^ \ {q + ie){k + ie){q +k +ie) 


xq+k+G;-{Q)G--{K). 

Again G~~{Q)G~~ (K) contain statistical functions = 
this by performing the following algebra: 


_ _ 4T^ _ 

qOfcO (2q~+q+){2k~+k+) 


4T2g+A:+ 


2^- 


_ rp2 _ 


T^q 


2u- 


T^k 


T^q-k- 


+ 


{2q'^+q-){2k++k~) 2g++g“ 2k++k~ (2g++g“)(2A:++A:“) 

The first term will not contribute: deforming the q'^ integral 


— adv 

(F.13) 
. We handle 

(F.14) 


5qi 


( 1 ) 


= g^GnGAT^ 


d^Q f d^K 


c^{2tt)^J {2tt)^ \{q -ie){k - ie){q +k - ie) 


1 + ^ 


2q+q- - 

(g+)2 ' q+J 2q+{q- - 5Eci + ie) 

the q~ integration can be closed below, yielding 
,4 


+ Ca ) 


dq 


+ ^ ^ Jcr (27r)3 J {2?:)^ \SEq{k--ie){k-+6Eci-ie) 


X 


(g+)2 


1 + 


q~ 


qI 


q+ J 2q+ 


+ Ca ■ 


Prr {CC) 

(F.15) 

Prr (-^) 

(F.16) 


The k integration can be closed in the upper half-plane, giving 


5qi 


( 1 ) 


= g^GnGAT^ 


dq^d‘^q± f dk^d‘^k± i 


X 


(<?+)' 


f dq+d^q^ f 

Jcr (2vr)3 J 

1 + QA+Ca. 


(27r)3 6El 


2q+ 


[Grik-=0)-G],-ik-=-6E^)] 


(F.17) 


This vanishes on Cr, because the square bracket is at least linear in dE^^. 

The second and third term are identical to Eq. (F.8) and thus vanish. Only the last 
term contributes, yielding 

44^ r ^4^ 27r%- + A:-) _ 


dqi 


= -g^GRGA 


f d^Q f 

J (27r)4y 


(2vr) 


-G--{Q)G--iK), 


(E.18) 


which cancels Eq. (E.ll). 
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Q + K 



Figure 14. The cat-eye diagram 


F.3 The cat eye diagram 

The diagram is shown in Fig. 14, and will be labeled c for cat-eye. Using only the bare 
vertex, the graph yields 


Sqi 


o^CrCa l^"dx+£ dx+' I ^ I 


xq'^{q~^ + k~^) 


G^%Q+K)G;,^{K)G-^^iQ) + G-/(Q+K)G/(K)G-^(Q) 


+G;P{Q+K)G-^iK)G-^{Q) 


(F.19) 


where we have defined the three-gluon vertex as 

gjabc-pf^up^p^ Q, K) = - QY + g’^PiQ - KY + g^^{K - PY] • (F.20) 

P,Q,K are all inflowing in the vertex, P is associated with a and /i and similarly for 
the others. The x"*" and x~^' integrals yield 


Sqr 




■J-. 


1 


{2ttYJ (27r)‘^\(g +i€){q +k Pie) 


- adv r^,p(-Q, -K, Q+K) 


xq^{q~^ + k^) 


G^YQ+K)G;YiK)GY^{Q) + gyYQ+k)GY'{k)Gy^{Q) 


+gyYQ+k)gyYk)g-^{q) 


(F.21) 


Let us look at the r/a structure of the propagators. Suppressing Lorentz indices and 
using Grr{K) = nB{K){GR{K) — Ga{K)), ng{K) ~ T/k^ , the last two lines can be 
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rewritten as 


Tq+Grr{K)[GR{Q)GR{Q+K) - Ga{Q)Ga{Q+K)] 
+T^p{Q)[Gr{K)Gr{Q+K) - GAiK)GA{Q+K)] 

-T^GrriQ) [p{K)GAiQ+K) + p{Q+K)Gr{K)] 

-r ^ GrriK) [p{Q)GAiQ+K) + p{Q+K)Gr{Q)] 

_TS_±^Grr{Q+K) [p{Q)Gr{K) - p{K)GRiQ)] 

+ ^ ^ S^GrriQ) [Grr {K)Ga {Q+K) + Gr {K)Grr (Q+K)] 

+Grr{K)GrriQ+K)GR{Q)]. (F.22) 


We start by dealing with the first line, which has the highest power of q~^ in the numer¬ 
ator. We label its contribution (1). It reads 


5qL 


( 1 ) 


= q^CrGa 


[ d^Q [ 

J (2vr)4 J 


d^K 


1 


(27r)^ V(g +iGj{q +k + ie 


— - adv) Tq+G-^{K) 
le) J 


y.V^,p{-Q,-K,Q^K) 


G/(Q)G/(Q+if) - G/(Q)G/(Q+i^) 


(F.23) 


Having obtained a fully retarded (advanced) function of q^ we can now expand on Cr 
{Ca)- Similar comments about pinching poles apply here as well: we expect Cf^{Q) and 
G"r{Q+K) to introduce poles for q~ = 5E^ and q~ + k~ = (5F^q+k respectively. Indeed 
we obtain 


dqi 


( 1 ) 


= q^GrGaT 


dq+ f dq-d^q± f d^K 


27r 


'Cr 

,2 , „ 1 , 


1 


(27r)3 J {2ttY \{q +*e)(g +k + ie) 

I 


— adv 


2{qi + qx • k^)GjiiQ)GjiiQ+K)G--iK) + O 


+ Ca- 


Expanding the transverse propagators to order 1 /( 5 +)^ we have 

GliQ) 


2q+{q -5Ec^ + iey 


Gl{Q+K) 


1 - 


q +k~ 


‘2q~^{^q -\-k — 6E^^]^-\-ie) y Q~^{q d~k — 6E^^]^-\-ie) 


(F.24) 


(F.25) 
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and considering only the (I/q’"*')® terms in Eq. (F.24) we have 


5qL 


( 1 ) 


= -q^CrCaT 


dq^ f dq d^q±_ f df'K 


'Cr 


27r 


(27r)" 


(27r)"^ \(q' +ie){q +k + ie) 


— adv 


(^1 + q± • {K) 


2{q+Y{q -(5Eq + ie)(g +k -(5Eq+k + ^e) 

+Ca- 


1 - 


q + k 


q+{q +k -(JEq+k+ie) 


(F.26) 


Rewriting the terms in round brackets as ^-functions gives 


dqi 


( 1 ) 


= -ig^CnCAT 


f dq~^ f dq d‘^q± f d'^K „ 1 /r/ - ,-x r/ -x\ 


(^i + q± • k L)G^^ {K) 


2{q+Y{q -5E^ + ie){q +k -(5Fq+k + *e) 
+C-A- 


where P denotes the principal value. This yields 


1 - 


q + k 


q~^{q Pk —5E^^\^+ie) 


(F.27) 


6qi 


( 1 ) 


= -ig^CRCAT 


xP-^ 


f f d‘^q± f d'^K {qj_ + qj_ • k_i_)G^r (K) 

Ic^ ^ J J 2{q+y 

1 k- 


+ - 


k [5Eq^{k -(iPq+k + ^e) 
+ Ca- 


1 - 


q+{k -dEq+k+ie) 


{k + JFq - ze)5Pq+k_ 

With some algebra and making pinches explicit we have 


(F.28) 


5qi 


( 1 ) 


= -ig^CRCAT 


X 


'Cr 

1 


dq'^ 


d^q^_ f d‘^K {q\ + q_L • k^)G^^ {K) 


f d^q± f 

J (2vr)2y 


(2vr) 


4{q+y 


+ - 


5Eci{k -5Ec^+k + ie) 

1 


k- 


' - A(k--5E,,^+.e)) [iTTTe + > 


—-^- 27ri6{k ) 


{k + 5Fq — ie)(5Fq+k \k — ie 
which gives, upon expanding the non-pinched denominators 


+ Ca, 


(F.29) 


5qi 


( 1 ) 


= -^g 


dq'^ f d‘^q± f d‘^K {q\ + q_L • k_L)G^r {K) 




X 


'Cr 

- 2TTi5{k~) + + 


(27r)4 2{q+f5E^5E, 


q+k 




{k~ + ieY {k~ — ieY 


+ Ca- 


(F.30) 


We have dropped terms that are 0(l/(g'^)^). The terms suppressed by one further power 
of l/c/"*" in Eq. (F.24) turn out to be either completely independent of q~ on Cr^ and 
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hence vanishing when its integration is done, or proportional to the terms in Eq. (F.26) 
times k~/q~^, q~jq^ or jq^. From the previous calculation it should be clear that 
the only way they could contribute at order 1/q'^ on Cr would be if both pinches {q~ 
and q~ + k~) were taken. In the first two cases that is not possible, because the factors 
of q~ or k~ at the numerator eliminate either of the two pinched poles and in the last 
case the resulting k'^ integration is odd once k~ is set to zero. 

We now consider the second line in Eq. (F.22), which we label (2): 


6qi 


xr^,p(-g, -K, Q+K) \G-j,^{K)G-^f’{Q+K) - G/(i^)G/(g+i^) 


We can now deform the k'^ integration, obtaining 


5qi 


( 2 ) 


= -q^GrCaT^ 


L 


dk~^dk d?k \ 


(2vr) 


d^g 

(2vr)4 


1 


, 2{kl + • ci^)Gl{K)Gl{Q+K)p 


{q +ie){q +k +ie) 

(Q) 


— adv 


/c+ 


(F.32) 


which is very similar to what we had before, due to the symmetries of the vertex. Higher- 
order terms in the expansion will not be relevant, as this contribution is one power 
smaller on the arc. Hence, replacing the transverse propagators with their leading-order 
expressions (F.25) and rewriting the terms in round brackets as d-functions we have 


Sqi 


( 2 ) 


+ig^GRGAT^ [ 
Jc 


Cr 

{kl 


dk~^dk d‘^k± 
k L • q±)/5' 


^2^P^(S(g- + k- 


■(Q) 


2{k+Y{k -5Ef^ + ie){k +q - dFlq+k + ' 


-Kd )) 

(F.33) 


which yields 


dqL 


( 2 ) 


p~~{-k-,q^,q±) _ 

(fe- - dFq + ie)dFq+k {k 


dq"^d‘^q± 1 -h k_L • q±) 

(27r)3 k- 2{k+f 

p (0,g+,g±) _ 

- (5Fq -I- ie){k- - dFq+k + 


(F.34) 


The second term on the bottom line vanishes under the q'^ integration, as it is odd. 
Similarly, the first term yields 


dqL 


( 2 ) 


-q^GrGaT^ [ 
JCh 


dk'^d?k± 

(27r)3 


dq'^d^q±_ {k\ + k_i_ • q_i_) 
(27r)3 2{k+)HE^+^5E^ 


X 


Gr (-(5Fq,g+,g_L) - (0,g+,g±) , 


(F.35) 
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which vanishes, as the integration can only pick up the residue of the Coulomb gauge 
poles, which is 0((5-Eq) and thus makes the k'^ integration vanish. 

Finally, terms with q~ or q~ + k~ at the numerator in Eq. (F.22) vanish again for 
the loss of q'^ at the numerator and of a pinched pole at the denominator. The last term 
trivially vanishes. The entire result is hence given by Eq. (F.30). 

F.4 Self-energy diagrams 

We analyze separately the two diagrams show in Fig. 15, the loop diagram on the left 
and the tadpole diagram on the right. 



Figure 15. The loop diagram on the left and the tadpole diagram on the right. 


F.4.1 The loop diagram 

The amplitude is labeled by I and reads 




P^CrCa f dq+d^q^ f d^K 


Tq+T^,p{Q, -Q - K, K)Tp,{-Q, Q+K, -K) 


2 J (27r)3 J (27r)4 

x<| + G'^/iQ+K)GPP'iK)] - adv 


q-=0 

(F.36) 


where 1/2 is a symmetry factor. We now perform a shift K ^ K — Q in the first term 
in square brackets on the second line. In principle one should be careful in performing 
such operations, as the integrals here are not finite. Indeed, as we anticipated, we will 
need to subtract the HTL counterterm, which, however, is obtained by performing the 
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same shift, as we will show in Sec. F.5. We then have 


Sqi 


X 


g^CnCA f dq+£q^ f d^K \ 
J (27r)3 J (27r)4 ^ \ 


. 2 J (2.)3 y 

r^,p{Q, -Q - K, K)rp,,,p, (-Q, Q+K, -K)G''j^'{Q+K)GP p'(K) 


(F.37) 


+rp,p{Q, -K, K - Q)T^,,,p,{-Q,K, -K + Q)G';i^'{K)GP p\k - Q) 


— adv 


q -=0 


We are now free to deform the contour, since in this case there are no statistical factor 
poles, as q~ is set to zero. We then have 

dq+ f d'^q± f d^K 2Tq+q\ 


dqL =9 ^CrGa 


'Cr 


xG^;^{K) - 


d4_ f d^q± r 

27r J (27r)2 J 

{k- + 2A:+) , J, 


(27r)4 iql+Ml)^[ 


{GliK + Q) + Gl{K-Q)) 


2q+ 


{Gl{K-Q)-G],iK + Q)) GY{K) 


+^- 


G5f(iF) /I (qx-kx)' 




2q] 




+ C-A- 


(F.38) 


5-=o 


Terms that had a linear term in the azimuthal angle at the numerator have been ne¬ 
glected. Since there are no pinching poles in k~ we can safely expand the K ziz Q 
propagators, yielding for their sum, up to order l/(g+)^ 


Gl{K + Q) + Gl{K-Q) = 


+ 


<^.E'q+k <JF'q_k 




TT6{k ) 

g+ ' 2q+ [(fc“+ie)2 ' {k- - ieY q+ " k-\' 

(F.39) 


whereas the difference is 

■yT 


G^(iF - Q) - Gl{K + Q) = + O 




(F.40) 


Plugging this back in Eq. (F.38) we obtain 

dq+ f d^q± f d^K 




5q, ^ 2g CrCaTj^^ j J ( 2^)4 (^2 + 


G7+W 


'ir6{k 


i 

+ 1^ 


'^-E'q+k ^.E'q-k 


+1- 


HTpT 

2\{k~+ieY ' (k~ — ieY g+ k~ 
.GVf{K) (I (qx-k^)2 


I +2k ) j. 1 


q^ 


2g^ k'^ 


+ 0(l/((?+)2)|+C^. 


(F.41) 


F.4.2 The tadpole diagram 

The amplitude, labeled by t, reads 

-ig*CRCA f dq+d‘^q± f d^K 


6qi 


fdq+^ f 

J (27r)3 J 


Tq^ 


2 J (27r)3 J (27r)4 

^ iy^qpugpa gppgva g^ag^p) adv 




(F.42) 


J<?-=o 
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Expanding on Cr we have 


dqi 


-ig^CRCAT f dq+ f (fq± f d^K 


-'Cr 


27r 




/c^ 


(27r)2 J (27r)^ 


G^j^{K) + 0 


q+{ql+Miy 


qlGYiK) 


(q+y- 


+ Ca- 


(F.43) 


F.4.3 Summary 

Summing Eqs. (F.41), (F.43) we obtain 


Sqi 


t+1 —ct 


= q^GrGaT 


f d^ [ 

Jcr 2vr J (27r)2 J 

B G- 


d^K 




+ 




{k~ — ieY 


i 


(27r)4 4.{q+f5El 
-{K) 


G7+(F')27r5(r 


{k~ — iey 


+ Gy^{K)[2- 


2k+k‘^. 


k'^{k~ — ie) 


-\-C.Ai 


(F.44) 


where we have used the fact that 5Ec^±\i = 5Fq + /c^/(2g+), up to vanishing terms in 
the azimuthal integration and that the square bracket on the second line, which can be 
identified with the NLO contribution to Zg in Coulomb gauge [27, 60], is purely real, so 
that the prescription used for the k~ poles at the denominator there is irrelevant. 


F.5 The subtraction term 


diflf 


We now turn to the computation of the subtraction counterterm dq^ , i.e. the soft 

subtr. 

part of the HTL self-energy. To this end, we need only the gluon loop, as soft fermions 
are not Bose enhanced and do not contribute to relative 0{g). The contribution from 
the loop diagram is, after the previously-discussed shift 


dqi 


difr 


1 subtr. 


g^GRGA f d^Q f d^K 


/ 


x<!G/(Q)G/ (Q) 


(2^)4 J ) 

{Q+K)Gi^J PP' {K) + {K)Gf {K-Q) 


— adv 
(F.45) 


where the vertices are treated in the HTL approximation, i.e. Q <C FT, and the 
propagators on the second line are bare and in the soft approximation, ng(k^) -A- T/k^. 
This yields 


dqL 


diff 
1 subtr. 


ig^GRGAT 


f dq+d^q^ 
Cr (2vr)3g+ 


rd^K gW^{k) 
y (27r)4 {ql+MlY 



(q±-kx)^ ' 

k^ 


+ Ca, (F.46) 


where we have used the fact that in Coulomb gauge the bare longitudinal spectral density 
vanishes and the transverse one puts K on shell. Furthermore, consistently with the HTL 
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approximation, in {K ± Q) one has only to keep the leading terms in K ^ Q, i.e. 

{Q ± K)'^ — >■ ±2Q ■ K, (q =t k)^ —)■ k'^. Other terms in the propagators and vertices do 
not contribute, as discussed in [27] . 

The tadpole contribution is trivially obtained from Eq. (F.43) by replacing the 
resummed K propagator with its bare counterpart, i.e. 

dq+d‘^q± f d^K r ■ -u .^2- 

so that the sum is 


Sqi 


diff 

= IO'^CrGaT I xg , 

tsubtr. Jr„ {2Tryq+ 


. f dq+d^q^ f 

L (27r)3g+ J 


(27r)4 (g2[' 


2 , (q± • kj 

1±+ k2 


+ Ca, (F.47) 


Sqi 


diff 


= ig^^GnGAT 

subtr. 


f dq+d^q^ f 

JCr (27r)3 J 


dq+d^q± f d^K 2q\Gf)'^{K) 


+ Ca- 


(F.48) 


F.6 


(27r)3 J (27r)4g+(g2 + 

Summary 

Summing the contributions from all diagrams and subtracting the counterterm (F.48) 
we obtain 

A.. ^ r rr d'K 
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loop 


- Sqi 
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+ - 
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G~~{K)TrS{k~ 
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where 5G = G - is the difference between resummed and bare propagators; in 
Coulomb gauge Eq. (F.48) is equivalent to the bare part of Eq. (F.44). After subtracting 
the collinear counterterm^® given in Eq. (5.12) and performing the K and integrations 
as in [27, 33] we obtain Eq. (6.8). 


G Semi-collinear integrations 

Let us consider Eq. (8.8). As mentioned, we put an IR cutoff on q_\_, which is the 

same cutoff used for the diffusion and conversion processes . We find it is simpler to use 
5E as an integration variable, so that, after performing the integration and the 5E 
integration with cutoff /{2\oj{p — a;)|) we have 
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shift in the integration variable is necessary, see footnotes 9-11 in [27]. 
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so that Eq. (8.8) turns into 
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Let us now explicitly obtain the logarithmic sensitivity to the diffusion region. Upon 
expanding Eq. (8.1) and its gluonic equivalent around cu ~ 0 we find 
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i.e. the expected diffusion structure (see for instance Eq. (3.23) or App. D) has appeared. 
Regulating the dix integral with an 0{T) UV regulator^*^ and expanding for small 
one obtains 
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which indeed cancels the p^^ dependence of Eq. (6.8). 

Similarly, as in Eq. (5.13), we can take the small-u; (or small p — oj) limit for the 
final-state quarks in Eq. (8.1) and examine the overlap with the conversion sector. We 
take as example the q ^ g rate, the opposite being the same times dp/dA- We obtain 
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Performing the same integration and expansion as before we then have 
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which removes the p^^^ dependence of Eq. (7.5). 


^Its exact value is not relevant, as we are only interested in the constant in front of the logarithm. 
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